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CHAPTER  1  -  INTRODUCTION 


A.  Background 

In  the  1960's  the  Navy  installed  the  first  mainframe  computers  to  manage 
their  vast  inventories  of  spare  and  repair  parts.  Along  with  these  computers  they 
installed  inventory  management  models  which  had  been  developed  by  Hadley 
and  Whitin  [3].  The  objective  function  of  these  models  was  the  minimization  of 
the  average  annual  total  variable  costs  to  procure  and  hold  inventories. 

The  Navy  manages  both  consumable  and  repairable  items.  Consumable 
items  are  discarded  when  they  cease  to  function  correctly.  For  repairable  items 
an  attempt  is  made  to  repair  a  nonfunctioning  item.  The  inventory  models 
developed  by  Hadley  and  Whitin  {3}  were  for  the  consumable  items.  Since  there 
was  no  model  in  Reference  [3]  for  repairables  the  Navy  decided  to  approach  the 
repairable  problem  by  subdividing  the  problem  into  two  distinct  parts,  those 
nonfunctioning  units  which  couldn’t  be  repaired  and  those  that  could  be 
repaired.  Each  part  was  "managed"  using  tiie  same  model  structure  as  was  being 
used  for  the  consumable  items.  The  units,  which  could  not  be  repaired,  were 
replaced  in  batches  through  a  procurement  action.  Those  that  coidd  be  repaired 
were  batch  inducted  for  repair. 

Using  the  two-part  approach,  the  Navy  was  able  to  develop  formulas  for 
the  economic  order  and  repair  quantities.  To  determine  the  two  reorder  points 
the  Navy  needed  to  have  a  backorder  cost.  However,  they  had  no  way  of 
determining  such  a  cost.  Therefore  they  adopted  an  approach  which  was  to  set  a 
goal  of  satisfying  an  average  annual  requisition  fiU  rate.  This  measure  was  called 
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the  "supply  material  availability"  or  SMA  and  the  overall  goal  was  an  SMA  of 
85%.  From  the  SMA  goal  an  implied  backorder  cost  could  be  determined. 

In  1982  the  Navy  attempted  to  integrate  the  two  parts.  Unfortunately,  the 
effort  was  only  partially  successful;  two  inventory  management  systems  still 
existed. 

In  the  late  1970's  the  decision  was  made  to  upgrade  the  mainframe 
computers.  The  Navy  decided  that  it  would  also  be  a  good  time  to  review  its 
models  and  improve  them  where  possible.  The  Naval  Postgraduate  School  was 
asked  to  participate  in  this  model  improvement  process.  In  1984  the  Navy 
accepted  a  wholesale  provisioning  model  by  Richards  and  McMasters  [7]  of  the 
Naval  Postgraduate  School  which  had  a  readiness-based  objective  function.  It 
was  the  minimization  of  the  aggregate  Mean  Supply  Response  Time  (MSRT)  for 
a  group  of  new  items  for  a  specific  weapon  system.  The  model  sought  maximum 
inventory  levels  for  the  items  which  would  satisfy  this  objective  subject  to  a 
provisioning  budget  constraint.  Unfortunately,  Navy  did  not  have  a 
replenishment  model  which  was  readiness-based.  Therefore,  the  provisioning 
model  was  never  used.  In  an  attempt  to  resolve  this  problem  tiiis  author  began 
the  search  for  a  replenishment  model  in  1986  which  has  the  same  objective 
function  as  ttie  provisioning  model  and  a  similar  budget  constraint. 

An  appropriate  model  for  managing  a  group  of  consumable  items  was 
developed  in  1989  [2].  A  model  for  managing  repairable  items  was  more  difficult 
because  of  the  complexity  of  the  process.  In  1988  a  study  of  preliminary 
simulation  resiilts  suggested  to  this  author  that  when  demand  is  modeled  as  a 
Poisson  process  that  the  probability  distribution  for  the  inventory  position  (on- 
hand  +  on-order  +  in  repair  -  backorders)  could  be  approximated  by  the 
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convolution  of  two  discrete  uniform  distributions,  one  for  repairable  carcasses 
and  the  other  for  carcasses  which  were  either  not  returned  or  not  repairable 
(attritions).  Batch  procurement  of  a  quantity  Qp  and  batch  induction  of  a  repair 
quantity  Qp  were  assximed  (see  References  [1]  and  [4]  for  the  details  of  the 
probability  distribution). 

The  simulation  model  required  furtiier  refining  based  on  discussions  willi 
operations  analysts  and  repairables  managers  working  for  the  Navy's  Inventory 
Control  Points.  The  current  form  of  the  simulation  model  for  repairables  was 
finalized  in  1992.  When  a  Poisson  demand  occurs  the  model  decides  if  there  is  a 
carcass  being  returned;  if  so,  tiien  that  carcass  enters  a  repair  queue;  if  not,  then 
that  information  is  sent  to  an  attrition  queue.  When  Qp  carcasses  have 
accumulated  in  the  repair  queue  the  entire  batch  sent  to  the  depot  for  repair. 
However,  they  are  usually  inducted  one  at  a  time.  As  each  is  inducted,  it  is 
examined  to  see  whether  it  is  capable  of  being  repaired;  if  not,  then  an  attrition  is 
added  to  the  attrition  queue.  Carcasses  which  can  be  successfully  repaired 
("good"  carcasses)  pass  through  the  repair  process.  The  first  "good"  carcass  goes 
immediately  into  repair  and  departs  a  repair  tumarotmd  time  (RTAT)  later.  The 
second  carcass  waits  a  short  period  of  time  (as  if  waiting  for  the  first  item  to 
finish  the  first  stage  of  repair)  and  then  enters  the  repair  process  if  it  can  be 
repaired;  otherwise,  it  is  rejected  and  recorded  as  an  attrition  and  the  next  carcass 
is  immediately  examined.  A  "good"  carcass  completes  repair  a  RTAT  later.  As 
each  "good"  carcass  in  repair  is  completed  it  is  returned  to  the  ready-for-issue 
(RFI)  inventoiy.  When  the  attrition  queue  reaches  a  size  Qp  a  procurement  of 
Qp  units  is  made  and  tiiat  batch  is  delivered  to  the  RFI  inventoiy  a  procurement 
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lead  time  (PCLT)  later.  A  flow  chart  of  this  model  is  presented  as  Figtire  1  in 
Chapter  6  of  this  report. 

Simvdation  results  from  this  model  [4]  showed  that  the  approximate 
inventory  position  distribution  developed  by  this  author  was  quite  robust;  that 
is,  it  accurately  represented  the  simulated  distribution  of  inventory  position  for 
any  set  of  Qp/Qj^  values  regardless  of  changes  in  values  of  all  the  other  model 
parameters. 

In  1994,  after  examining  the  simulation  results  of  the  repairable  item 
inventoiy  management  process  just  described.  Baker  [1]  was  able  to  apply 
stochastic  modeling  techniques  to  derive  an  approximate  formula  for  the 
probability  distribution  for  the  value  of  the  net  inventory  at  any  instant  of  time. 
This  allowed  the  next  step  in  the  model’s  development,  which  was  to  derive  the 
formulas  for  tiie  determining  of  the  probability  of  being  out  of  stock  at  any 
instant  of  time  and  the  expected  number  of  backorders  at  any  instant  of  time  for 
a  repairable  item.  The  importance  of  knowing  the  expected  number  of 
backorders  at  any  instant  of  time  is  that  it  can  be  used  to  determine  the  MSRT. 
The  probability  of  being  out  of  stock  at  any  instant  of  time  is  needed  to 
determine  the  expected  number  of  backorders  at  any  instant  of  time  and  the 
SMA. 

Two  other  aspects  of  the  model  have  also  been  studied;  safety  stock  and 
economic  order  and  repair  quantities.  The  Navy  uses  the  safety  level  (safety 
stock  quantity)  to  provide  a  measure  of  protection  provided  by  an  inventory  of 
an  item.  For  a  repairable  item  there  exists  no  theoretical  formula,  although  the 
Navy  does  have  a  procedure  for  determining  one  since  they  need  it  for 
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developing  reorder  points  for  the  two  parts  of  their  model.  A  simulation  study 
was  conducted  this  past  year  to  determine  an  approximate  formula  for  safety 
stock  which  would  apply  for  the  process  modeled.  The  definition  of  this  safety 
stock  [4]  is  that  it  is  "the  expected  net  inventory  at  the  time  an  order  of  new  imits 
arrives  into  the  RFI  inventory  and/ or  a  repaired  carcass  is  returned  to  tiie  RFI 
inventory."  Four  approximate  formulas  were  fotmd  which  fit  the  simulated 
results  quite  well;  two  for  the  case  of  batch  repair  (no  delay  between  "good" 
carcasses  being  inducted),  and  two  for  the  case  of  "good"  carcasses  being 
inducted  one  at  a  time  with  a  delay  between  inductions.  The  results  of  that 
study  are  presented  in  this  report. 

The  formulas  for  the  Navy’s  economic  order  and  repair  quantities  are 
modifications  of  the  dassic  square-root  formula  known  as  the  "economic  order 
quantity"  which  is  designed  for  use  with  consumable  items  [3].  They  do  not 
change  when  the  reorder  point  is  changed  although  there  is  an  iterative 
procedure  in  Reference  [3},  Chapter  4  for  doing  so.  The  question  which  came  to 
mind  when  thinking  about  the  development  of  a  readiness-based  model  was 
what  effect  does  a  changing  maximum  inventory  position  have  on  the  economic 
values  of  the  order  and  repair  quantities.  Do  they  change  or  remain  fixed?  A 
study  was  therefore  conducted  to  answer  ttiis  question.  The  results  of  that  study 
are  presented  in  this  report 

B.  Objectives  and  Scope 

The  first  objective  of  this  report  is  to  present  the  derivations  of  the 
formulas  for  determining  the  probability  of  being  out  of  stock  at  any  instant  of 
time  and  the  expected  nrimber  of  backorders  at  any  instant  of  time  for  a 


5 


repairable  item  for  the  net  inventory  probability  distribution  derived  by  Baker. 

In  addition,  based  on  discussions  with  personnel  at  the  Navy’s  Inventory  Control 
Point  (NAVICP),  a  modification  of  the  assumption  of  Poisson  demand  during  the 
aggregate  lead  time  to  a  Normally  distributed  demand  was  requested. 

Therefore,  the  second  objective  is  to  present  the  derivations  of  the  formiUas  for 
that  modification.  A  third  objective  is  to  present  the  results  of  a  simulation 
analysis  conducted  to  discover  approximate  formulas  for  safety  stock.  Finally, 
the  fourth  objective  is  to  present  tiie  results  of  economic  analyses  conducted  to 
determine  approximate  formulas  for  the  economic  order  and  repair  quantities. 

C.  Preview 

Chapter  2  presents  the  derivations  of  formulas  for  the  probability  of  being 
out  of  stock  at  any  instant  of  time  under  the  assumptions  of  Baker’s  model. 
Chapter  3  presents  the  formulas  for  the  expected  number  of  backorders  at  any 
instant  of  time  tmder  the  assumptions  of  Baker’s  model.  Chapter  4  presents  the 
derivations  of  the  formulas  for  the  probability  of  being  out  of  stock  at  any  instant 
of  time  imder  the  assumption  of  a  Normal  demand  during  lead  time.  Chapter  5 
presents  the  derivations  for  the  formxalas  of  the  expected  number  of  backorders 
at  any  instant  of  time  under  the  assxunption  of  a  Normal  demand  during  lead 
time.  Chapter  6  presents  the  details  of  the  simulation  study  used  to  derive 
approximate  formulas  for  safety  stock  for  repairable  items.  Chapter  7  presents 
the  details  of  the  development  of  approximate  formulas  for  determining  the 
order  and  repair  quantities  which  minimize  the  average  annual  total  variable 
costs  of  managing  an  inventory  for  any  specified  value  of  the  maximum 
inventory  position.  Chapter  8  presents  a  brief  summary  of  the  previous  chapters. 
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conclusions  with  respect  to  the  chapter  results,  and  reconunendations  for  further 
model  development 
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CHAPTER  2  -  DERIVATION  OF  PouriSW)  FOR  BAKER'S  MODEL 


A.  Introduction 

This  chapter  presents  the  derivations  of  the  formulas  for  the  probability  of 
being  out  of  stock  at  any  instant  of  time  based  on  the  net  inventory  distribution 
developed  by  Baker  [1].  To  facilitate  the  process,  Baker’s  distribution  is  first 
converted  to  a  form  involving  complementary  distribution  functions  of  the 
Poisson  distribution.  This  conversion  then  allows  the  use  of  identities  from 
Appendix  3  of  Hadley  and  Whitin  [3]  in  the  derivations  of  the  PQjjj-iSW) 
formulas. 

B.  Conversion  of  the  Probability  Mass  Function  of  Baker's  Model 

Equation  (38)  from  Bako:'  s  thesis  [1]  for  the  case  where  demand  during 
the  aggregate  lead  time  is  Poisson  was 
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QpQrPo 

min(Qp,Qj^)  ^ 
QpQk  ,-."i 

1  minrar^,zK 
I  ^  y  <^Mfly 

QpQr  j=x^+\ 


+ 1  - 
(z-j)! 


JbrO  ^2  <Xi 


forxi<z  <X2 


forx2  <  z; 


(2.1) 


where 


x-^  =  miniQp,Qp)-\ 


^2  ^Max  (2.2) 

^Max  ~  Qp  ■*■  Or  “  2 ' 

N(t)  is  the  net  inventory  at  any  instant  of  time,  and  SW  is  the  maximum 
inventory  position.  As  shown  in  Baker's  thesis,  fi  =  Z  where  Z  is  the  Program 
Problem  Variable  or  the  modified  form  ZB  which  indudes  REP.  The  format  of 
equation  (2.1)  is  difficult  to  use  to  derive  formulas  for  Pqijj(SW )and  the  time- 
weighted  xinits  short,  B( SW) . 

The  following  form  facilitates  the  use  of  identities  from  Appendix  B  of 
Hadley  and  Whitin  [3].  In  this  form  P(.)  is  the  Poisson  complementary 
cumulative  distribution  function  and  p  is  the  mean  of  the  distribution. 
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kz  +  l-n)-(z  +  l)P(x>z  +  l)  +  fiP(x>z)  for  0<z<Xi; 


(xj  +l)-(z  +  l)P(x>z  +  l)  +  nP(x  >z) 

+(z  -Xi)P(x  >z-Xi)-^iP(x>z-x-^-l)  forxi<z<X2; 


QpQrP(SW-z) 


-(z  +  l)P(x  >z  +  l)  +  fiP(x  >  z) 
■\-(z-Xi)P(x>Z-Xi)-fJLP(x  >Z-Xi  -1) 
+(2  -X2)P(x>Z~X2)-fjP(x  ^Z-X2-V 
-^Max  -  ^  2  -  inin(Xj^^  ,Z)) 

+p,P(x  >z- min(xj^x /Z) “ D 


forx2<z. 


(2.3) 


The  derivation  of  equation  (2.3)  follows.  Beginning  with  the  first  term  of 
equaticm  (2.1) 


S 


(z-j)l 


/ 


letx  =  z-|;  then)  =  z  -  x  andj  +1  =z-x  +1.  For  ttie  limits  of  the  summation, 
when  j  =  0  then  x  =  z  and  when  j  =  z  then  x  =  0.  Asa  consequence. 


!AL_Ll=y(z_^+i)iL£ 
po  ~ 


x! 

z 


=  (z+l)^p(x;fi)-'^  xp(x; n), 
x=0  a:=0 

where  p(x;  p)  is  the  Poisson  probability  mass  function;  that  is. 


p(x;p)  = 


x! 


(2.4) 


Notice  that  the  summation  limits  when  x  replaces  z  -  j  are  reversed.  The 
summation  upper  limit  for  j  corresponds  to  the  lower  limit  for  x  and  visa  versa. 
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This  equivalence  can  be  seen  by  expanding  the  two  summations.  When  this  is 
done  we  get 


I 


/=o 


ife  ^  2n^  V  ^  ^  (z-l)juV^  z^e~^  (2  +  l)e“^ 

2!  (2-1)!  (2-2)1  2!  1'  0!  ' 


and 


£(2-I  +  1)L 

x=0 


(2  +  l)e  ^  2//g  ^  (2-l)jU^e  ^ 

0!  1!  2! 


3/i 


+  -T - -  +  - - 


(2-2)1  (2-I)! 


2I 


The  first  term  of  the  second  line  of  equation  (2.4)  can  be  easily  rewritten  as 

(2  +  1)X  /^)  =  (2  +  ^P{x  >  0)  -  P(x  >  2 + 1] 

x=0 

=  (2  +  l)[l.0-P(x>2  +  l)], 
since  P(x  >  0)  =  1.0  for  the  Poisson  distribution. 

The  second  term  of  the  second  line  of  equation  (2.4)  can  be  simplified 
using  Identity  1  from  Appendix  3  of  Hadley  and  Whitin  [3] 

2  Z  2 

X  xp(x;  p)  =  X  =  p^p(x-l;n)  =/x[1.0  -P(x>z)]. 

a:=0  ;c=0  2=0 

Finally,  then, 

2  /  •  ,  Z—j  -/i 

(  J  =(2  +  l-/l)-(2+l)P(x>2  +  l)  +  /lP(y>2).  (2.5) 

y=o  ~ 
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For  the  range  of  <  z  <X2  the  probability  mass  function  given  by 
equation  (2.1)  adds  a  term  and  changes  one  of  the  boimds  on  the  summation  of 
the  first  term.  The  first  term  is  now 


^  =(2+1)  £  xpix;^) 

7=0  '  ^  ’  X=2~Xi  X=Z—Xi 

=  (z  +  l)[P(x  >z-x{)-P{x^z+D]  (2.6) 

- p,P(x  >  Z  -Xj  -1)  +  pP(x  >  z) . 

The  second  term  is 


miniQp,QR)  X  . 

;=2i+l 


=  (^1+1)  S 

:e=0 

=  (xi  +  l)[l.0-P(x>z-Xi)]. 


(2.7) 


Note  that  min(Qp,Qj^)  - 1  =  Xj  (see  equation(s)(2.2)) 

Summing  equations  (2.6)  and  (2.7)  and  collecting  like  terms  yields 


=  (xj  + 1)  -  (z + l)P(x  >  z  + 1)  +  juP(x  >  z)  (2.8) 

+  (z-Xi)P(x>Z-  Xi)-^P(X>Z-Xi  -1). 


For  z  >X2,  the  first  line  is  the  same  as  for  Xj  <  z  ^  X2.  The  second  line  has 
the  upper  limit  as  Xj  instead  of  z. 


mm(Qp,Qg)  £  + 


>=^1+1 


(2-y)! 


x! 


(2.9) 


=  (Xi+l)[P(x>Z-X2)-P(x>Z-Xi)]. 
The  third  line  takes  a  littie  more  work.  First, 
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nun(2,ar^) , 


^ (^Mox  +1  - 7)/^^  ^  iXMax+l-Z+x)n^e  ^ 

^^2+1  x=z-Tmn(z,x^  ) 


2-^-1  Z-Xj-1 

=  (^Mflx  +  l-2)  X  Pix;tl)+  ^  xpix;^). 

x^z-  Tmn(z,XMax )  j:=2-min(2,XM« ) 

Applying  Identity  1  of  Appendix  3  of  Reference  [3]  leads  to 
i^Max  +1  - 


mn{z,x^) , 

I 

;=x,+l 


(2-7)1 

=  (^Max  + 1  -  z)[Pix  >z-  min(z,XMax)  -Pix>z-  X2)] 

Z-X2  -1 

+  n  ^  p{x-l;n) 

x=z-inin(z,Xj^) 

=  i^Max  +'i--z)[P(x>  Z-  inin(2,Xjv^) -P{x>z- X2)] 


(2.10) 


+  /f[P(x  > 2  -min(2, X^ax)  -1) -P{x>Z-X2-  1)] • 

Summing  equations(2.6),  (2.9),  and  (2.10)  and  collecting  like  terms  gives 
the  density  function  for  2  >^2  in  the  revised  form  shown  in  equation(2.3). 

y=o 


2  —  +  »™(ep^Qf!)  L  ‘t—KT*  E 


/=xi-n(2-;)l  '  ;=t'+l 


=  -(2  +  l)P(x  ^  2  + 1)  +  llP{x  >  z) 

+  (z-Xi)P(x>Z  -  Xi)  -  flP(x  >Z  -  Xi  -1) 

+  (2-X2)P(x>2-X2)-7iP(x>2-X2  -1) 

-  (z-  Xmox  -1)P(x>z-  minCxM^^, 2))  +  fiP(x  <  2 -  min(xA^„^ ,2) - 1) . 


C.  Derivation  of  PoujiSW) 

1.  The  case  of  0  <  SW  <  Xi 
The  general  formula  for  Pout(^^) 
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(2.11) 


Pom(SW)=  £  [N(t)  =  SW-z)]. 

z=SW 

Note  that  when  z  =  SW  then  N(t)  =  0 ;  that  is,  there  is  no  stock  on  hand. 
Obviously  then,  there  will  be  no  stock  on  hand  when  z  >  SW ;  the  system  will  be 
in  backorder  status  with  the  number  of  backorders  being  given  by  z  -  SW . 

The  next  step  is  to  apply  tihe  equation  (2.3)  results  to  equation(2.11).  We 
begin  with  the  case  for  0  <  S W  <  Zj. 


Poiir(SW)  = 


1 

QpQr 


Xi 

-pL)-(z  +  l)P(x>Z  +  l)+  liP(x> 

2=SW 


z)] 


+ 


1  ^  (a:i+l)-(z  +  l)Px>z  +  l)+/iP(x>z) 

QpQr  2=^+1  -Xi)P(x  >  z  -Xi)-nP(x  >z  -Xi  -1) 


-(z  +  l)P(x  >z  +  l)  +  ^(x>z) 


+ 


1 

QpQr 


oo 


E 


+(z-X^)P(x>Z-Xi)-ljP(x  >Z-Xi-1) 
+(z  -X2)P(x>Z-X2)-fiP(x  >Z-X2-1) 
-(x -Xmox  - '^)P(x  ^z- min{z,xj^ }) 
+HPix  >  z  -  min{z,XMax  }  - 1) 


(2.12) 

This  can  be  further  simplified  by  collecting  common  terms.  For  example, 

-(z  +  l)P(x  >z+1)+illP(x>z) 

appears  in  all  three  parts  of  equation  (2.12).  Therefore,  for  this  case  it  can  be 
summed  for  z  from  SW  to  <» ;  tihat  is, 

oo 

£  [-(z  +  l)P(x  >  Z  +  1)  +  flP(x  >  z)]. 

z=SW 
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Next,  the  tenn 


(z-Xi)P(x>z-Xi)-^(x>z-Xi-l) 

shows  up  in  the  second  and  third  parts,  so  it  can  be  summed  for  z  from  arj  to  oo. 

oo 

[(z-Xi)P(x>z-Xi)-^iP(x>z-x^-l)]. 

z=jei  +1 


The  last  terms  come  from  the  third  part  for  which  the  upper  summation 
boimd  for  z  is  already  oo.  From  titiere  we  have 


(Z-X2)P(X>Z-X2)-^IP(X  >Z-X2-1) 


1 


-(z  -^Max  -VP(x  >  2  - min{z,xj^^}) 
+^(x  >2  -  tnin{z,Xn^ax}  -1) 


There  are  also  the  constant  terms.  From  the  first  part  we  have 

i  i  2+  i  a-n) 


z=SW 


z=SW  z=SW 

x^(xt  +  1)  (SW-l)SW' 
2  2 


+(1-  ^i)[xi-  max(0,SW  -1)]. 


For  the  second  part  we  have 

(X-^+l)=(Xi+l)(X2-Xi)  . 

2=Xj  +1 

Thus,  we  can  rewrite  Pout 
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Poiir(SW)  = 


QpQr 
+ 


x^(x-^+l)  (SW-l)SW 
2  2 


1 

Q^^[(l-fi^Xi  -  max(0,SW  -1)]  +  (x^  +  l){x2  -Xj)] 
1 

7777“  S  [-(z  +  DPfx  >z  +  lJ  +  /iPfA:>z)] 

WpWk  z=SW 


2,  [(Z  -^1  )P(^  ^  2  ^  z  -Zj  - 1>] 


QpQr  2=r+l 


-J—  T 

QpQi?  z=t+i 


(Z-X2)P(X>Z-X2)-HP(X  >Z-X2-1) 

-(2  -Xmox  -  VP(x  ^  z  -  min{z,XMax  }) 
[+tiP(x  > z  -  mm{z,ZA^^}  - 1) 


(2.13) 

The  first  two  lines  of  equation  (2.13)  is  consists  of  the  constant  terms 
(i.e.,  no  probabilities  need  be  evaluated).  The  rest  of  the  equation  will  require 
identities  from  Reference  [3]  to  simplify  the  form  of  the  summations. 

The  simplest  form  is  the  last  term  in  the  third  line;  namely. 


fiP(x>z)  =  ti  2,  P(x>z)=^a(SW). 


2=SW 


z=SW 


(2.14) 


For  convenience,  we  have  used  the  following  shorthand  notation  in 
equation  (2.14); 

oo 

a(u)=  '^P(x>z), 
z-u 

and,  from  Identity  6  of  Appendix  3  of  Reference  [3],  it  follows  that 

a(u) =iiP(x>u-\)-(u- \)P(x  > u) .  (2.15) 
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As  a  consequence  of  the  definition  of  a(u)  it  follows  that  the  other  terms 
in  equation  (2.13)  involving  can  be  readily  written  as  forms  of  aiu). 

oo  oo 

-  2  ^lP{x>z-Xl-l)  =  -^l  X  P{x>z-x^-l)}  (2.16) 

Z=Xi  +1  Z=Xj+l 

Let  y  =  2 -ACj-l.  Then,  if  2  =  ^1 +l,y=0.  If  2  =  «>,y  =  «5.  Therefore,  we  can 
write  equation  (2.16)  as 

oo  oo 

-H  X  P{x>z-Xi-l)  =  -nJ^Pix>y)  =  -tia{0). 

2=;Ci+l  y=0 

Likewise,  we  can  let  y  =  2  -  22  - 1  in  the  last  part  of  the  fifth  line  and  get 

00  00 

-PL  2  P{x>z-  X2-1i)=-p^P(,x^y)  =  -iLoiO).  (2.17) 

z=X2'¥\  y=0 

From  the  definition  of  a{u)  it  follows  that 

a(0)  =  iiP{x  >-l)  +  P{x  >  0)  =  JU+ 1, 
since  P(x  >  -1)  =  P(x  >  0) = 1.0  for  the  Poisson  distribution. 

The  next  term  to  be  simplified  is  the  first  term  in  the  third  line  of  equation 
(2.13).  The  process  first  sets  y  =  2 + 1 .  Then,  when  2  =  SW,  y  =  S  W  + 1 . 
Therefore, 

00  00 

2  (z  +  l)P(x>2+l)=  2^  yP(x>y)  =  y(SW  +  l), 

z=SW  y=SW+l 

where,  using  Identity  8  of  Appendix  3  of  Reference  [3], 

r(«)=  ^yn.xty)  =  L^P(.x2u-2)*iiHx>u-l)-!^^^Pi;ci:u).  (2.18) 

y=«  ^  ^ 

When  we  let  y  =  2  -  Xj  and  y^z-x^,  respectively,  equation  (2.18)  allows  us  to 
write 
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oo  OO 

X  (z-A:i)P(x>z-Xi)  =  Xy^(^-y)  =  /(!)' 

z=x,+l  y=l 

OO  oo 

X  (z-A:2)PU>2-X2)=Xy^(^-y)  =  )'(l)*  (2.19) 

Z=X2+1  T/=l 


Two  terms  still  remain  to  be  simplified.  These  are  ttie  sbcih  and  seventh 
lines;  namely, 

oo 

X  ~^Max  -  ^)P(x  >  z  -min{z, Xmox}) ; 

Z=X2+1 


X  tiPix^z-mm{z,XMax}-'0- 

2=a:2+1 


It  is  useful  to  separate  the  summation  into  two  parts;  from  ^2  + 1  to  Xj^,  and 
from  xj^  +1  to  oo.  The  first  summations  result  in 

'£i 

z=X2+l  z=X2+l  2=a:j+l 

~  (^Max  ■*"  ~  ^^Max  1)^1  / 


(^Max  +  'i-)  2-  H^^2-mm{z,Zjviax})=UAto+l)  Z  P(x  >  0)  +1) 


X  HP{x>Z-mm{z,XMax}-'^)  =  ^^  S  Pix>-1)=  n{XMax-X2)=  m  ) 

2=a^+l  z=X2+l 


and 
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-  ^  zP(x>z-min{z,XMax})  =  -  S  zP(x^ 


0) 


z=x,+l 


Z=Xj+l 


^Miu  r 

=  -  I  — 

2=^2 +1 


^Max  (^Max  +  ^2  ^^2  + 


]■ 


Putting  these  parts  together  as  they  are  in  equation  (2.13)  gives 

^Max 

-  X  zP(x>z-min{z,XMax}) 

2=Zj+l 

2=Zj+l 

+  X  ^2-wm{z,x^^}-l) 


2=2^  +1 
“^l("^5Wa2  +  1  +  f^)  ~ 


^Max (^Max  1)  ^2(^1  "t  1) 

2  2 


The  final  part  is  to  simplify  these  components  for  the  summations  from 
^Max  ■*■1  to  00  . 


00 

X  (-Z  +  Xmox  +1)P(x>z-  min{z,XMax}) 

2=2j^+l 

00  00 

—  ”  ^  ^  ^  ^  Z  —  ^Max^  (2.20) 

OO  OO 

=  -'^}/P(x'^y)+'^P(x>j/)  =  -Y(l)  +  a(l). 
y=i  y=i 


Here  we  set  y  =  z-Xf^a^.  When  2  =  x^^  +  lthen  y  =  l. 


20 


Next, 


5^  flP(x>Z-mm{z,XMax}-V  =  ^  £  P(^^^-^Max-V 

2=Xm^+1 


Z=X*to  +  l 


(2.21) 


=  /i'^P(x>y)  =  fia(0). 
y=o 

For  this  part  we  set  y  =  z-xj^^ -Iso  when  z  =  + 1  tiien  y=0  . 


Finally,  combining  these  two  parts  gives 


-  5^  zP(x>z-min{z,XMax})  +  (xMax  +  V  S  P(x>z-min{z,XMax}) 

Z=^Mu+l 


2=Z*to+l 


(2.22) 


+  5)  /iP(x^z-mm{z,Xjvifl;t}-l>=-y('l)  +  a(l)  +  /^a(0). 


z=z«„+l 


In  conclusion,  the  sum  of  the  last  two  lines  of  equation  (2.13)  is 

oo  oo 

-  2  (z-XMax-VP(x>z-min{z,Xj^ax})+  X  MP(x^z-mm{z,XMax}-l) 

2=X2+1  2=:i:2+1 

=  Xi(xMax+'^+M)-  j -yd)  +«(!) + Ata(0; . 


Earlier  we  showed  a  simplified  form  for  a(  0) .  However,  for 
programming  purposes  it  is  convenient  to  leave  it  in  the  o!(«  jgeneral  format. 
That  will  also  be  true  for  y(l) .  The  definitions  for  ofw)  and  y( u)  were  given  in 
equations  (2.15)  and  (2.18),  respectively. 
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QpQr 


x^(x^+l)  (SW-1)SW 
2  2 

H-Cl  -  ^l)[xl  -  max(SW  - 1)]  +  (x^  +  l)(x2  -x^ ) 
-r(SW  +  1)+  na(SW)  +  yd)  -  na(0) 


+a(\) 


^Maxi^hAax  +  V  _^2fa2  + 
2  2 


(2.23) 


Note  that  summing  the  resulting  simplifications  of  equations  (2.17),  (2.19),  and 
(2.20)  gives 

-na(0)  +7(1)  -  7(1)  +a(l)  +/ia(0)  =  a(l) 
for  the  last  line  of  equation  (2.23). 

2.  The  case  of  Xi  <  SW  <  X2 

The  next  case  for  Pqjjj(SW)}  namely,  that  for  Xj  <  SW  <  X2,  has  only  the 
second  and  third  parts  of  equation  (2.13)  and  ttie  lower  boimd  on  the  second  part 
is  now  SW. 


WSW)  = 


1 

QpQr 


“I 

s 

2=SW 


(Xi  +l)-(z  +  l)Px>z  +  l)  +  ^(x  >z) 
+(z-Xi)P(x>Z-Xi)-^(x  >Z-Xi  -1) 


-{z+l)P(x  >2  +  1)+  fjP(x>z) 


+ 


1 


+(2-21  )P(X  >  2  -Xj  )  -  flP(x  >  2  -Xj  - 1) 
+(Z-X2)P(X^Z-X2)-^(X  ^Z-X2-V 
“(Z  -XMax -'i-)P(x  >  2  - min{z/)CMax}) 


.  (2.24) 


+fiP(x  >  2  -  min{z^^gx  }  - 1) 
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The  constant  tenn  in  the  first  line  is 


X  fXi  +  D=('xi+lXr2-SW+l).  (2.25) 

z=SW 

Again,  for  the  second  and  third  components  of  the  first  line  we  have 

oo 

X  [-(z  +  l)P(x  >z+l)+ nP(x  > z)]  =  -y(SW +1)  +  na(SW).  (2.26) 
2=SW 

since  these  parts  are  also  the  fliird  line. 

The  second  line  of  equation  (2.24)  is  now 

oo 

^  [(z-Xi)P(x>z-Xi)-^iP(x>z-Xi-l)] 

2=SW 

y=SW-x^  y=SW-Xi-l 

=  y(SW  ~Xi)-iia(SW -Xi-l). 

Equation  (2.27)  includes  the  third  line  of  the  second  part  of  equation  (2.24). 

The  last  three  lines  of  equation  (2.24)  do  not  change  from  the  final  form  of 
Pout(^^ )  for  ttie  case  where  0  <  SW  <  that  is. 


+(z-X2)P(x>Z-X2)-flP(x  >Z-X2-1) 
-(z -Xmox  - VP(x >z- mm{z^Max}) 


^  [+jUP(x  >  2  -  min{z,Xj^  }  - 1) 

= 7  (1)  -  ^la(0)+Xl(X}^ax  +  l  +  ^l)- 


XMox^X^ijx  1)  _^2(^2 
2  2 


-7(l)+a(l)  +  jua(0) 

XMox^Xmox  ***  22(^2  1) 

2  2 


=Xi(XMax+'^  +  P-)- 


+a(l). 
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Putting  all  the  pieces  together  gives,  for  Xj  + 1  <SW  <  X2, 


Pqiji(SW)  = 


QpQr 


(Xi  +  1)(X2  -SW  + 1)  -  y(SW  +  1)  +  iia(SW) 
+y(SW  -Xi)-^a(SW  -x^  -1) 

+ 1  +  ji)-'  >‘2(^2  + V 

+a(l) 


.  (2.28) 


3.  The  case  of  X2  <  SW  <  Xj^^gx 

In  this  case  all  that  remains  of  equation  (2.3)  is  the  third  part. 


Pout(^^)- 


QpQr  z=SW 


-(z  +  l)P(x  >Z  +  1)+  nP(x  >  z) 
+(z-Xi)P(x>Z-X-^)-ljP(x  >Z-Xi  -1) 
+(z-X2)P(x>Z-X2)-UP(x  >Z-X2-1) 
-(^  ~^Max  - '^)P(x  ^Z- min{z,Xj^a^ }) 
+HP(x  >z-  fnin{z^Max}  - 1) 


.  (2.29) 


There  are  no  constant  terms  in  equation  (2.29).  The  summation  results  for 
the  first  two  lines  have  been  determined  above  in  equations  (2.26)  and  (2.27).  AH 
that  remains  to  be  derived  are  the  summations  for  the  last  three  lines  of  equation 
(2.29).  Similarly  to  the  derivations  leading  to  equations  (2.17)  and  (2.19),  we  have 


00  00 

-/i  X  P(x>z-X2-1)  =  -H  X  P(x>y)  =  -^(x(SW- X2-I) ; 

2=sw  y=SW-Xj-l 
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^  (z-X2)P(x>Z-X2)=  X  }fP(x>}/)=Y(SW- X2). 

z=SW  y-SW-X2 


Two  terms  still  remain  to  be  simplified.  These  are  the  last  two  lines; 
namely, 

00 

-  X  (^-XMax-'i-)P(x^z-min{z,XMax})! 

z=-SW 

00 

X  ^P(x>z-  min{z,XMax} "  1)  • 

z=SW 


It  is  useful  to  separate  the  summation  into  two  parts;  from  SW  to  x^,  and 
from  Xjviox  +1  to  00.  The  first  summations  result  in 

(XMax+'^)  X  P(^^^-^H^'^Max})=(XMax+'^)  X 

z=SW  z==SW 

“  (^Max  ^  / 


X  IlP(x>z- 

2=SW 


min{z,xj^}-l)  =  fi  X  P(x>-l)  =  ^i(xMax-SW+l); 
z=SW 


and 


-  X  zP(x>z-min{z,XM^})=-  X  zP(x>0)  =  -  X  2 
2=SW  z=SW  2=sw 

^  \xMax(XMax  +  V  _(SW-1)SW 

I  2  2 
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Putting  these  parts  together  gives 


-  ^  zP(x>z-min{z,Xj^^„^}) 


z=SW 

+(^Max  +  S 

2=SW 


+  S  ^lP<x 


z=SW 


>z-mm{z,XMax}-V 


(^Max  1  "I"  t^)(^Max  + 1) 


XMax(XMax  +  l)  (SW -1)SW 
2  2 


Equation  (2.22)  provides  the  restilts  for  the  summations  from  Xj^  +1 
to  oo.  Thus,  for  the  case  where  X2  <  SW  <  x^^  the  sum  of  the  last  two  lines  of 
equation  (2.29)  is 


oo 

-  Y,  (z-XMax-'l^)P(x^z-mm{z,XM^}) 
z=SW 

oo 

+  ^  ^lP(x >z- min{z, Xj^^x } - D 
2=SW 

=  +1  +  ^^)(XMax  -  SW  +  1)  -  _(SW^^)SW 

-Y(V+a(l)  +  fia(O). 
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Therefore,  when  X2  <  SW  < 


PoiiT(sy^)= 


1 

QpQr 


_  ^Max(^Max  + 1)  fSW  -1)SW 
[2  2 
^-r(SW  +  l)  +  ija(SW) 
+r(SW-x-^)-fia(SW-XT^  -1) 
+y(SW  -X2)-fia(SW -X2  -1) 

-Y(l)-i-a(l)  +  ^ux(O) 


(2.30) 


5.  The  case  of  SW  > 

For  the  case  where  SW  >  x^^j^  the  first  two  lines  of  equation  (2.30)  are  no 
longer  applicable  and  the  last  two  lines  of  equation  (2.29)  are 


X  (-z  +  XMax+l)P(x>z-min{z,XMax}) 

z=SW 

00  00 

=  -X  (z-XMax)P(x^Z-XMax)+  X  P(x^Z-XMax) 

z=SW  z=SW 

00  00 

=  -  X  yP(x>y)+  X  P(x>y)=-r(SW-XMax)+a(SW-XMax); 
y=SW-x^  y=SW-x^ 


and 


00  00 

X  jUP(x>Z-mm{z,Xjvia;c}-l)=M  X  -1) 

z=SW  2=51^ 

00 

=  /r  X  P(x>y)=iia(SW-x^ca-^)- 

y=SW-x^-l 
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Therefore,  for  SW  > 


PqIjj(SW)  = 


1 

QpQr 


'-y(SW  +  l)  +  fia(SW) 

+Y(SW -Xi)-fia(SW -Xi  -1} 
'  -^yiSW  -Xi)-  na(SW  -X2  -1) 
-y(SW  -XM„^)+a(SW  -Xj^) 
+fux(SW —1) 


(231) 


D.  Summaiy 

This  completes  the  derivations  of  the  equations  for  Pqijj(SW ).  Equation 
(2.23)  applies  when  0  <SW  <  Xy  equation  (2.28)  applies  for  x^  <  SW  <  Xi, 
equation  (2.30)  applies  when  X2  <  SW  <  x^ax'  equation  (2.31)  applies  for 

SW  >  Xj^ax  •  These  are  the  equations  which  were  programmed  for  the  Poisson 
distribution  representing  the  demand  during  aggregate  lead  time  for  a  repairable 
item. 
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CHAPTER  3  -  DERIVATION  OF  B(SW)  FOR  BAKER'S  MODEL 


A.  Introduction 

This  chapter  presents  the  derivations  of  the  formulas  for  all  of  the  cases 
for  the  expected  number  of  backorders  at  any  instant  of  time,  J5(SW),  for  the 
Baker  model  [1]  where  demand  during  aggregate  lead  time  is  Poisson.  The 
derivations  will  rely  heavily  on  the  derivations  in  Chapter  2.  Indeed,  we  begin 
with  the  case  where  0  ^  SW  <  Xj  and  recall  equation  (2.12).  Instead  of  merely 
summing  all  elements  of  the  net  inventory  probability  distribution  provided  by 
Baker  (equation  (2.2),  which  the  modified  form)  as  we  did  for  Pqijj{SW),  we 
must  include  the  number  of  backorders  for  each  case;  that  is,  z-  SW  when 
2  >  S  W .  However,  we  will  also  include  the  case  where  z  =  S  W  so  that  we  can  use 
the  derivations  of  Pqijj{SW)  as  part  of  the  backorder  derivations.  We  do  realize 
that  when  z  =  S  W  that  there  are  no  backorders  but  it  is  mathematical  convenient 
to  include  it  anyway.  For  Baker’s  model  we  define  B(SW)  as 

oo 

B(SW)=  X  (z-SW)[N(f)  =  SW-2].  (3.1) 

z=SW 

Next,  we  can  separate  this  equation  into  the  following  two  parts; 

oo  oo 

B(SW)=  X  z[N(f)  =  SW-2]-SW  X  [m)  =  SW-z]  .  (3.2) 

z=SW  z=SW 

The  second  term  is  nothing  more  than  the  product  SW Pquj(SW).  Therefore  we 
can  focus  ovu:  attention  on  the  derivations  associated  with  flie  first  part.  We  will 
denote  this  part  as  Bi(SW). 
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(33) 


Bi(SW)  = 


I 


z=SW 


iN{t)  =  SW-z]. 


Therefore,  we  can  write 


BiSW)  =  BiiSW) -SWPouriSW) .  (3.4) 


B.  Derivation  of 

1.  The  case  of  0  <  SW  <  Xi 

For  the  case  of  0  <  S  W  ^  Xi,  when  the  net  inventory  probability 
distribution  is  given  by  equation  (2.2)  we  can  write  Bi(SW)  as 


Bi(SW)= 


1 

QpQr 


1 


2=SW 


+  1-  ^)-iz+l)P{x>z  +  l)+  ^iP{x>z)] 


+ 


1 

QpQr 


ixi+l)-{z  +  l)Px>z  +  l)  +  fiP{x>z) 
+iz  -  Xi)Pix  >Z-Xi)-tjP{x>Z-Xi-l) 


(3.5) 


-(z+  l)P(x  >  2  +  1)  +  nP{x  >  z) 


^  1  y  ^ 

QpQRz=x^+1 


+(z-Xi)P(x>Z-Xi)- flP(x>Z-Xi  -1) 
+(Z -X2)P(X>Z-  X2)  - /lP(x  > 2  -  X2  - 1) 
-(2  -  Xmox  -  Wx  >  2  -  Tc6xi{z,XMax  }) 


+/xP(2  >  2  -  mm{2,ACjvto}- 1) 


Next,  following  the  form  of  equation  (2.12),  equation  (3.5)  can  be  rewritten  as 
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Bi(SW)  = 


QpQr 


2  z(z+l-)U)+  ^  z(Xi  +  l) 

Z=SIV  Z=Xi  +  l 


1 

+  7r7r  S  2[-(2+l)B(:s:>z  +  l)  +  /iP(a:>z)] 

2=SIV 

+  77^'  S  z[(z-Xi)P(x>z-a:i)-/xP(x>z-Xi-1)] 

WJPWR  2=Xj+l 

{Z-X2)P{X>  Z -X2)- HPix>Z- X2-I) 

TTTT  2  -(z -  X^ax  -  Wx  >Z- iran{z,XMax }) 
+liP(x  >  2  -  mm{z,2M«}- 1) 


(3.6) 


Addressing  the  constant  terms  first, 

^  z(z  +  l-/i)=  Z^+(1-M)  i  2 


2=5^ 


z^SW  z=SW 

Zi(a:i  +  l)(2yi  + 1)  SW(SW  -1)(2SW  -1) 
6  6 
1(1  ^)'^’(^'  +  l)  (SW-1)SW 


(3.7) 


.Z  JC 

2  2(Zi  +  l)  =  (Xi+l)  Z=(Xi  +  l) 


Z=Xi+l 


Z=Zi+l 


X2{X2+^i  zi(yi+l) 

2  2 


.  (3.8) 


The  second  line  of  equation  (3.7)  is  a  consequence  of  the  identity 

^  ^2  _  w(w  +  l)(2w  + 1) 

2=0  ^ 


Adding  equations  (3.7)  and  (3.8)  and  collecting  like  terms  gives: 
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2(z+l-/i)+  ^  2(2^  +  1)  = 

2=SW  2=r  +  l  ■  L  ^ 


+1)  SW(SW-1)(2SW-1) 


i  ■ 


^1(^1 +  1)' 


•M 


2i(2t  +  1)  SW(SW-1)~ 
2  2 
2;  (22+1)  5W(SW-1) 


(3.9) 


We  next  need  to  develop  tiie  parts  involving  the  probabilities.  The  third 
line  of  equation  (3.6)  is  first. 


^  2[-(2  +  l)P(x>2+l)+/fP(x>2)] 


2=SW 


(3.10) 


=  -  ^  Z^P{x>Z  +  l)-  ^  zPix>Z+l)  +  tl  ^  zP(x>z). 
2=SW  2=SIV  2=SIV 

Looking  at  die  first  term  on  the  right-hand  side,  let  y  =  2  + 1 .  Then, 

z2=(y-l)2=y2-2y  +  l. 

Also  when  2  =  SW  ttien  y  =  S  W  + 1.  This  change  of  variable  gives 


-  ^  2^P(x>2  +  l)  =  -  ^  (/ -  2y  +  l)P(x  >  y) 

2=sw  y=SW+l 

=  -  E  y'^Hx>y)  +  2  X  yP(x>y)-  X  P^x>y). 

y=SW+l  y=SlV+l  y=SW+l 


(3.11) 


Next  we  define  5{u)  =  X  3^  ^(^  -  y^  Identity  9  of  Appendix  3  of  Reference{3} 

y=u 

provides  the  following  result: 
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5(u)='f^/Pix>y)=  P{x>u)+tiPix>u-l) 

^  -*  (3.12) 

2  3 

+  ^-P{x>u-2)+^P{x>u-3). 

The  definitions  of  aiu)  and  /(u)  were  given  in  Qxapter  2  by  equations 
(2.13)  and  (2.16),  respectively.  Thus,  the  second  line  of  equation  (3.11)  in  final 
reduced  form  is: 

oo 

-X  z^P{x>z  +  l)  =  -diSW  +  D+2r{SW  +  l)-a{SW  +  l).  (3.13) 

2=SW 

The  second  term  of  equation  (3.10)  is  similarly  analyzed.  We  get 
zP(x>z+l)=-  X  {}f-l)P{x>y) 

z=m  y=SW+l 

=  -  X  }/Pix>y)+  X  P(x>SW)  (3.14) 

y=SW+l  y=SW+l 

=  -y(SW+l)+a(SW+l). 

The  third  term  of  equation  (3.10)  is  next.  No  transformation  of  variable  is 
needed. 

oo 

H  X  zP(x  >z)-  jayiSW) .  (3.15) 

2=SW 

The  final  form  for  equation  (3.10)  is  the  sum  of  the  right-hand  sides  of 
equations  (3.11),  (3.12),  and  (3.13): 
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X  z[-(2  +  l)P{x  >  z+ 1)  +  ^iPix  >  z)]  =  -S{SW+ 1)  +  YiSW  + 1)  +  uriSW).  (3.16) 
z=SW 


Next  we  look  at  the  third  line  of  equation  (3.6).  It  can  be  subdivided  into 
three  parts. 

oo 

X  2[(2-  Xi)P{x  ^z-Xi)-  nPix  ^Z-Xi-1)] 

Z-X^-hl 


=  X  z^Pix>z-Xj) 

z=z,+l 

oo 

-Xj  X  zPix>z-Xi) 

z=x^+l 

oo 

-fl  X  zPix>z-Xi-l). 

Z=Zi+l 

The  first  term  on  the  right-hand  side  can  be  reduced  as  follows: 


(3.17) 


X  z^Pix>z-x{)  =  ^iy+x^fPix>y) 

z=z,+l  y=l 

OO  oo  oo 

=  X  ^  y)  +  2xi  X  1/P(x  >y)+ Xj^X  ^  y)  (3.18) 

y=l  y=l  y=l 


=  5(1) + 2xiy(l)  +  Xi^a(l) . 

The  second  term  is 

oo  oo 

-Xj  X  zP(x>z-Xi)  =  -XjJ^(y+x^)P(x>y)  =  -x^Y(D-Xi^a(l).  (3.19) 

z=Zi+l  y=l 

The  third  term  is 


-fi  X  2P(x>2-Xi-l)=-ju£(y  +  Xi+l)P(x>y) 

Z=a:i+1  y=0  (3.20) 

=  -HY(0)-n(xj  +  l)a(0). 


34 


Summing  equations  (3.18),  (3.19),  and  (3.20)  results  in  equation  (3.17) 
becoming 

oo 

y  Z[(Z -  Xi)P(x  >  Z  -  Zj)  -  llP{x  >  Z  -  Zj  - 1)1 
2=tVi  (3.21) 

=  5(1) + Ziy(l)  -  iiyiQi)  -  ju(zi  +l)a(0) . 


The  next  term  to  consider  is  the  first  line  of  the  last  part  of  equation  (3.6). 

oo 

y  z[(z  - X-^P{x  >  Z  -  X2)  -  jUP(x  >  Z  -  JC2  - 1)] 


(3.22) 


=  X  Z^P{x>Z-X2) 

2=^2+! 

00 

-X2  y  ZP(X>Z-X2) 
2=22+1 
00 

-)U  X  zP{x>Z-X2-l). 
z=x,+l 


Following  the  steps  for  equation  (3.18)  we  can  write  the  first  term  on  the  right- 
hand  side  of  equation  (3.22)  as 

OO  00 

y  Z^P(x>Z-X2)='^(j/  +  X2fP{X>y) 

2=2j+l  y=l 

=  '^/P{x>y)+ 2X2X  yP(x  >y)+  X2^  y  P{x  >  y)  (3.23) 

y=l  y=l  y=l 

=  5(1) + 2x27(1)  +  X2^a(l) . 


Following  the  steps  of  equation  (3.19)  we  can  write  the  second  right-hand  term  of 
equation  (3.22)  as 
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-X2  X  zPix>z-X2)=-X2'^{y+X2)Pix>y)  =  -X2r(l)-X2Wl).  (3.24) 

z=X2+l  y=l 


Finally,  following  equation  (3.19)  we  can  write  the  third  right-hand  term  of 
equation  (3.22)  as 


-H  X  zP(x>z-X2-l)=-n'^{y+X2  +  DPix>y) 

2=Xj+l  y=0 

=  -HYi0)-^i{x2  +  l)a{0), 


(3.25) 


Summing  the  right-hand  sides  of  equations  (3.23),  (3.24),  and  (3.25)  results 
in  equation  (3.22)  being  reduced  to 


X  Z[(Z  - X2)Pix  >Z-X2)-  llPix  >  Z  - X2  -  1)] 


(3.26) 


=  5(1)+X2y(l)-  /ry(0)-  /r(A:2  +l)a(0) . 


The  last  two  lines  of  equation  (3.6)  are 

oo  oo 

X  -z[iz-XMax-'^)Pix>z-Tr^iz,Xj^„^)]  +  ^i  X  zP(x>z-mm{z,XMax}-'0 

Z=X2+1  z=x^+l 

As  in  the  PQUjiSW)  derivations,  we  subdivide  this  summation  into  two  parts; 
from  X2+1  to  Xj^ax  ^Max  +1  fo  In  the  first  part  of  the  stunmation 

interval  we  have 


X  -z[{z  -  Xf^ax  - 'OPix  >  2  -  rcM.Z,X}^ax)] 


2=:c,+l 


=  X  -Z^Pix>0)+{XMax+'^)  X  zPix>0), 

Z=X2+1  z=x.  +  l 
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since  z - Min{z, }=z-z  =  0.  Likewise, 

=  /x  zP{x>-l). 

2=^:2+!  z=X2+l 

Because  P(ar  >  0)  =  1.0  and  P{x  >  -1)  =  1.0  for  the  Poisson  distribution  these 
equations  are  reduced  to 


fl  X  zPix>z-mm.{z,XMax}-^) 


2  -z[{z -  X^ax  - ^)Pix  >z- min(2,a:^^)] 

z=x^+l 

^Ata 

=  X  -2  +(^M«C  +  1)  X  2 

z=X2+l  2=:x:2+1 


(3.27) 


and 


+  (^Max  +M 


^Max  ^  Max 

fl  2  zPix>z-ir^{z,Xj^a^}-l)  =  fi  2  2 
z=X2+l  z=Xj+l 

(3.28) 


Combining  equations  (3.27)  and  (3.28), 


£  -z[(z  -  X^ax  -  VP(x  >  2  -  })] 

2=^2 +1 

+  M  X  zP(x^z-mm{z,Xj^gx}-l) 

2=^2+i  (3.29) 

__  ^Maxi^Maz  +  ^)(^Max  +1)  ^2(^2+ 1)(^2  + 

6  6 


+  (^Ato  +1  +  /^) 


^Max(^May  ^2(^2 

2  2 


The  remaining  part  of  the  summation  is  for  x  + 1  to  oo . 


£  -z[{z -  Xmox  -l)Pix>z-  min{z, x^ax })] 

2=X*^+l 

OO 

+  ^^  X  ^  “ min{z, - 1) 

2=**to+l 

OO 

=  £  -z[(z  -  Xi^^  -  1)P(2:  >  z  -  Xmox] 

oo 

■*■  M  X^  XP(X  >  Z  —  —  1)  . 

2=^ai«+1 


(3.30) 


The  first  line  of  the  right-hand  side  of  this  equation  can  be  reduced  as 
follows: 
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X  -4^^-  ^Max  - ^  2 -  XMax  ] 

2=a:*ta+l 

oo  OO 

—  ~  X  2  P(Ar  >  Z— Xj^^)  + +  1)  ^  zP{x>Z  — 

oo  OO 

=  “X  +  ^Max  fP(x  >  y)  +  + 1)5^  (y  +  Xmox  Wix  >y) 

y=i  y=i 

= -X  y  -  y ) "  X  y  ^  y )  -  X  ^  y ) 

y=l  y=l  y=l 

oo  oo 

— y)  — y) 

y=l  y=l 

—  —5(1)  —  ^^Max  l)y(l)  1)^(1) 

—  —5(1)  —  —  l)y(l)  +  ^Max^(^)  • 


The  last  line  of  equation  (3.30)  is  reduced  to 

oo  oo 

H  2;  zPix>Z-XMax-'i)  =  tl'^(}/+XMax+^)Pix>y) 


(3.31) 


x=x^+l 


y=o 

=  fiyio)  +n(xMax  +  l)a(0) . 


(3.32) 


This  completes  the  derivations  for  Bi(SW)  for  the  case  of  0  <  SW  <  Xj. 


Equation  (3.6)  now  has  the  following  form: 


39 


Bi(SW)  = 


QpQr 


xi  (xi  +  lX2x^  + 1)  SW(SW  -  1)(2SW  - 1)' 
6  6 
X2(X2+1)  XiUt+I)' 

2  2 


+a: 


x-iix-i+l)  SW(SW-l)' 


2  2  ] 

X2(X2  +  1)  5W(SW-1) 

"^2  2 

-5(SW  + 1)  +  y(SW + 1) + ^lYiSW) 

+5(1)  +  a:i7(l)  -  ny{0)  -n{xi  +  l)a(O) 

+5(1)  +  ^27(1)  -  uriO)  -  ^i{x2  +  l)a(0) 


XMaxi^Mtix  +  '^)i'^XMax  +1)  X  2(^2  +  1)(2X2  + 1) 
6  6 

XMaxiXMax+'^)  X2(X2+D 
2  2 


-5(1)  -  (^Alax  -  l)r(l)  +  XMaxfxil) 


+M7(0)  +  +  l)a(O) 


.  (3.33) 


Note  that  there  are  terms  which  can  be  cancelled;  such  as  5(1)' s  in  the 
seventh  and  tenth  lines  and  the  ///(O)  terms  in  the  seventh  and  eleventh  Unes. 
All  of  the  7(1)  terms  can  be  collected  and  reduced. 

^l7(l) X^ (1)  ~(^Ato  ~  1)7(1)  ~  [^1  ~  (^Alax  “^2)  ■tl]7(l) 

= [^1-^1+ i]r(i)= 7(1)- 

Similarly,  the  a(0)  terms  can  be  collected. 
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-fiix-i  +I)a{0)-^{X2  +  I)ai0)  + fi{x^  +l)a(0) 

=  [-Xi  -  3^2  +  Xmox  - 1  - 1  +  l]jt^a(0) 

=  [-Xi  +  Xi  -l]n(x{0)  =  -fia{0) . 

This  allows  consolidation  of  the  sixtiv  seventh,  tenth  and  eleventh  lines  into 


S{l)+XiYO)-nriO)-nixi  +l)a(0) 

+  5(1) + X2r{l)-HYiO)-^iix2  +  l)a(O) 
~  5(1)  —  (^jvto  ~  l)y(l)  + 

+  juy(0)  +  M3:Mflx+l)a(0) 

=  5(1) + 7 (1) +XAtoa(l)  -  MO)  -  ;i«(0) . 


Equation  (3.33)  is  then  reduced  to 


Bi(SW)= 


QpQr 


+x 


-At 


xr(x^  +  l){2x^  +1)  SWjSW  -1)(2SW  -1)' 

6  6 

'^2  (3^2  +  1)  3:i(Xi+l) 

2  2 

x^jx^+l)  SW(SW-l) 

2  2 

JC2(a:2  +  l)  SlV(SW-l) 

"^2  2 

-5(SW  + 1)  +  7(S  W  + 1) + HYiSW) 

r^Mfl3:(^Mgj:  ^)(2'^Mgj;  ^)  ^2(^2  '^l)(2r2  ~^^) 

1  6  6  . 

^Mg3:(^Mfly  ~^^)  ^2  (^2 't'  1) 

2  2 


■K^Mgx  +1+ At) 


[+5(1)  +  7(1)  +  XAto«(l)  -  At7(0)  -  HaiO) 


(3.34) 


However,  from  a  programming  point  of  view  little  is  gained  by  using 
equation  (3.34).  In  addition,  as  the  forms  of  Bi(SW)  for  the  other  cases  are 
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developed  below  it  is  useful  to  compare  them  witii  equation  (3.33)  rather  than 
equation  (3.34). 


2.  The  case  of  Xi  <  SW  <  X2- 


Equation  (3.6)  is  now  reduced  to 

Bi(SW)=-^  i 

WPWR  2=:SW 

■‘■TTTT  S  2[-(2+l)P(a^^2  +  l)  +  /fP(ac>2)] 


^  z[{z-x^)P{x>z-Xi)-nP(x>z-x-i-l)] 


QpQr  z=SW 


QpQs  ^-7*1 


iz-X2)P(x>Z~X2)- flPix>Z-X2  -1) 
-(z-  Xmox  -  '0P{x  >z-mm{z,XMax}) 
+HP{x  ^  z  -  mm{z,Xf^ax}- 1) 


(3.35) 


The  constant  term  is  easily  reduced. 

z=sw  z=m  ^  ^  ^ 


.(3.36) 


The  second  and  fourth  through  sixth  lines  of  equation  (3.35)  remain  the 
same  for  this  case  as  for  0  <  SW  <  x^.  The  third  line  dianges.  Following  the 
equation  (3.17)  analysis,  we  have 
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z[(z - x-i)P(x  >  z-  Zj) - flP{x  >  z -Xj -  1)] 

2=SW 

=  X  z^Pix^z-Xi) 

z=SW 

oo 

-^1  2  zP{x>z-Xi) 
z=SW 

oo 

zP{x>z-Xi-l). 

z=SW 


(3.37) 


The  first  term  of  equation  (3.37)  is  reduced  as  follows; 

X  2^P(x  >z-Zi)=  (y  +  Xi)^P(z>y) 

2=sw  y=SW-2, 

OO  oo  oo 

=  2  /P(^^3/)+2xi  5^  yP(x>y)+xf  ^  P(x>y)  (3.38) 

y=SW-Xi  y=SW-a;i  y=SW-x^ 

=  5(SW-  xi)  +2:ciy(SW  -Zi)+zfa(SW  -xj) . 


The  second  term  becomes: 


oo  oo 

-Xj  2  2P(2>z-Xi)  =  -Xi  X  {y+x-^)P{x>y) 

z=SW  y=SW-x,  (3.39) 

=  -Xiy(SW  -  Xj)  -  xfa(SW  ~x^). 


And,  the  third  term  is: 

oo  oo 

zP(x>z-x^-l)=-/i  £  (y+Xi+l)P(x>y) 

2=sw  y=SW-a:,-l  (3.40) 

=  -/j.  y(SW  -  X|  - 1)  -  fi(xi  +  l)ce(S  W  -  x^  - 1) . 
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Summing  the  results  of  equations  (3.38),  (3.39),  and  (3.40)  gives  equation 
(3.37)  in  the  following  reduced  form: 

oo 

z[iz- Xi)Pix  >z-Xi)- nP{x  >  z  - - 1)] 

z=SW 

=  5(SW-Xi)+Xiy(SW-Xi)  (3.41) 

-  fiyiSW-x^  - 1)  -  nix^  +l)a(SW -x-^-l). 


The  form  for  Bi(SW)  is  therefore  given  by  equation  (3.42)  for  this  case. 


Bi(SW)= 


(Xi  + 1) 


X2(X2  +  1)  (SW-1)SW 
2  2 


QpQr 


-6(SW  + 1)  +  y(SW  + 1)  +  hy(SW) 

+5(SW-Xi)+Xir(SW-x^) 

-fir(SW-x-^  -l)-H(x^+l)a(SW-x^  -1) 

+5(l)+X2Y(l)-fiY(0)-fi(x2  +  l)a(0) 

\^Max(^Max  +  ^X^Max  +  VC^X^  +  1) 

6  6 


+(^Max  +  1  + 


^Max(^Max 

2  ~  2 


- (^Max  -  l>>r 


+juy(0)  +  piUxi^ax  +  l-to:(0) 


(3.42) 
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Examination  of  equation  (3.42)  shows  that  the  5(1)  and  juy(0)  terms  will 


cancel.  In  addition,  the  y(l)  and  a(0)  terms  can  be  combined: 

-/X(X2  +  l)a(0) + }x{xMax  +  l)a(0)  =  -  ^2)«(0)  =  /i^ia(O); 

X2y{i)-{xj^  -i)r(i)=[i-(^Ate  -^2)]y(i) =(1-3^1  )r(i)- 
Equation  (3.42)  can  then  be  somewhat  reduced  to 


Bi(SW)  = 


QpQr 


(Xi+1) 


X2{X2  +  ti  (SW-1)SW 
2  2 


-5(SW  + 1)  +  y(SW  +1) + iiyiSW) 


+5(S  W  -  Xj)  +  Xiy{SW  -  Xi) 


-fiyiSW  -  Xj  - 1)  -  n{xi  +  l)a(SW  -  Xj  - 1) 
Ml  -  3Ci)7(l)  +  fix^aiO)  +  Xj^axO^O) 


^Maxi^Max  ^^2(^2+ 1)^^2  + 1) 

6  6 


M.xm„  + 1+ 

.  1_  z  z  ^ 


(3.43) 


2.  The  case  of  X2  <  S  W  <  x^ax 

In  this  case  we  need  to  consider  only  the  third  part  of  equation  (3.5). 
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Now  BiiSW)  has  the  form: 


Bi(SW)  = 


QpQr  z=SW 


-{z+l)P{x  > Z  + 1)  +  nP(x >  z) 
-^{z-Xi)P{x>Z-Xi)-^(x>Z-X^-l) 
■Kz-X2)Pix  >Z- X2)-pP(x  >Z-X2-1) 

-<z- Xmox  -1)P{x>z- min{z, x^ax }) 
+HP(x  >  2  -  minfz,  x^^x  }  “  D 


.(3.44) 


The  constant  terms  from  the  top  lines  of  equations  (3.33)  and  (3.42)  are  no 
longer  present.  The  final  form  of  die  first  line  of  equation  (3.44)  was  provided  by 
equation  (3.16).  The  final  form  of  the  second  line  was  given  by  equation  (3.41). 

The  derivation  steps  for  equation  (3.41)  provide  a  short  cut  for  the 
derivation  of  the  final  form  of  the  third  line  of  equation  (3.44).  When  X2  replaces 
Xi  in  equation  (3.38)  we  get 

X  z^P{x>z-X2)=  X  {y+X2fP{x>y) 

z=SW  y=SW-x^ 

=  X  y^P{x>y)+2x2  X  X 

y=SW-Zj  y=SW-X2  y=SW-a:j 

=  5(SW~  X2)+2x2r(SW-X2)+xia(SW-X2) . 


Similarly,  equation  (3.39)  jdelds 

00  00 

-Xz  X  zP(x>z-X2)==-X2  X  (y+X2)P(x>y) 

z=SW  y=SW-x^ 

=  -XzYiSW  -Xz)-  3^aiSW  -Xz). 
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And,  finally,  equation  (3.40)  gives 


zPix>z-X2-l)^-fi  2)  (}f  +  X2  +  l)P{x>y) 

z=SW  y=SW-a:2-l 

=  -fiyiSW  -  a:2  - 1)  -  Ii{x2  +  l)a(S  W  -  ^2  - 1) . 


Then,  combining  these  results, 

oo 

2)  2[(z -X2)P{x  >Z-X2)- llP{x  >  2 - ^2  - 1)] 

2=51^ 

=  5(SW-  X2)+X2y{SW -X2)  (3.45) 

-  /X7(SW  -  X2  - 1)  -  ix{x2  +  l)a(SW  -  JC2  - 1) . 


The  last  two  lines  of  equation  (3.44)  need  to  be  analyzed  for  the  two 
intervals  2:2 + 1  ^  S  W  <  X}^  and  S  W  >  Xj^  + 1.  The  last  two  lines  for  the  first 
interval  are  split  into  two  summations;  those  fi-om  SW  to  and  from  X}^,^  +1 
to  00.  FromSW  to  Xj^ax' 


^  -z[{z-XMax-\)P{x>Z-VC^{z,XMax)] 

z^SW 


+  fl  2  zP(x>z-mm{z,XMax}-'0 

z=SW 

=  X  -Z^Pix>0)  +  {XMax  +  T^)  L  zP{x>0)+n  £  zPix>-D  (3.46) 


2=sw  z=SW  2=SW 

'XMariXMar  +l)(2XMar  +1)  SW(SW-  1)(2SW-  1) 
6  6 

+  (XMax  +1  +  /i) 


47 


Equations  (3.31)  and  (3.32)  provide  the  remainder  of  the  summation  from 


^Max  +1  to  oo.  Therefore,  when  X2  <  SW  <  x^^x'  Bi(SW)  is 


B^(SW)  = 


-^(SW  + 1) + y(SW  +  1)  +  hy(SW) 


+6(SW -Xi)+Xiy(SW  -Xi)- ^iy(SW  -Xi  -1) 


-/f(jC|  +  l)Df (SW  -Xi  -1) 


+5(SW  -X2) + X2Y(SW  -3C2)  -  juy  (SW  -^2  -1) 
1  I  -ii(x2  +  l)a(SW  -X2  -1) 


QpQr 


’xMax(XMax  +  V(^Max  +  V  SW(SW  -1)(2SW  -1) 
6  6 

^Max(^Max  + 1)  S  W(S  W  - 1) 

2  2 


+ 1  +  A^) 


-d(l)-ixj^ax -1)y(V+X}^{1) 


+juy  (0)  +  ii(xj^ax  + 


.  (3.47) 


4.  The  case  of  SW  >  Xm^^ 

When  SW  >  Xj^^,  Bi(SW)  is  simplified  even  further  because  the  constant 
terms  of  equation  (3.47)  which  are  the  fourth  and  fifth  lines  are  eliminated.  The 
last  line  now  involves  S  W .  As  a  consequence,  the  lower  boimd  on  the 
summation  in  equation  (3.31)  is  now  SW  instead  of  +1  as  shown  in 
equation  (3.48) 
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“5(SW  —  “  ^Max)  ""  “  ^Max) 

(^Max  +l)y(SW  +  ^Max^^Max  l)Of(SW-^  ^Max) 

=  —5(5 W  —  Xj^^  )  —  {Xj^fnx  l)y  (S W  —  ^JVfox  )  (S W  —  Xj^Q^  )  . 
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Equation  (3.32)  provides  the  structure  for  the  second  line  of  equation(3.48). 


H  5^  zPix>z-XMax-l) 
z=SW 

oo 

=  H  X  (y+XMax  +  'i)Pix>y) 

y=SW-x^-l 

=  HriSW  -  X^ax  - 1)  +  tlixMax  +  l)a(S W  -  Xj^^  - 


The  final  form  for  Bi(SW)  when  SW  >  Xj^^x  ^ 


B^(SW)  = 


QpQr 


-d(SW  + 1)  +  y(SW  +  1)  +  nr(SW) 

+5 (SW -Xj ) + Xiy(SW  -Xi)- /iy (SW -x^  - 1) 
-H(xi+l)a(SW-Xi-l) 

+d(SW  -X2)  +  X2Y(SW  -X2)-tlY(SW -X2-V 
-ju(x2  +  l)a(SW  -X2  -1) 

S(SW  -X^ax)  - (xmox  -  Vy(SW  -Xj^^ ) 
’^Xj^ojficiSW  ~Xj^i^ )  +  ^Y  (SW  —  1) 

^P^X^fax  +  l)cr(SW  ~Xj^^  —1) 


(3.50) 


(3.51) 
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C.  Summary 

This  completes  the  derivations  for  Bi(SW).  In  summary,  for  B^CSW), 
when  0  <  SW  <  Xi  we  use  equation  (3.33)  or  (3.34);  when  Xj  <  SW  <  X2  we  use 
equation  (3.42)  or  (3.43);  when  X2  <  SW  <  Xj^  we  use  equation  (3.47);  when 
SW  >  Xf^g^  we  use  equation  (3.51).  Finally,  to  get  B{SW)  we  use  equation  (3.4) 
and  the  corresponding  formulas  for  Pquj(SW)  from  Qiapter  2. 
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CHAPTER  4  -PouriSW)  DERIVATION  FOR  THE  NORMAL 

MODIFICATION 

A.  Introduction 

The  derivations  presented  in  this  chapter  foUow  those  of  Chapter  2. 
However,  the  assumptions  of  discreteness  for  Inventory  Position  and  Net 
Inventory  are  no  longer  used.  As  a  consequence,  we  will  first  present  the 
formulas  for  the  Inventory  Position  and  Net  Inventory. 


B.  Inventory  Position  Distribution 


The  first  step  is  to  redefine  the  probability  distribution  for  the  Inventory 
Position.  The  Inventory  Position  distribution  for  the  discrete  demand  model  was 
the  convolution  of  two  discrete  uniform  distributions  (see  Baker  [1],  pp.  16, 17). 
Then,  using  the  logic  of  Hadley  and  Whitin  [3]  in  their  Chapter  4  (page  192)  we 
now  have  a  distribution  which  is  the  convolution  of  two  continuous  uniform 
distributions;  namely. 


/(^p)  =  ^  forO^Xp<Qp; 


f(xR)=-^  JbrO<Xp<Qp} 


(4.1) 


where  xp  is  the  number  of  attritions  in  the  procurement  batching  process  at  time 
t  and  Xp  is  the  number  of  carcasses  in  the  repair  batching  process  at  time  t. 


Letting  x  =:Xp+Xp,it  follows  that  if  SW  is  the  maximxun  Inventory  Position  and 
I(t)  is  the  Inventory  Position  at  time  t,  then  I(t)  =  SW  -  x . 

The  distribution  of  I(t)  is  then 
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X 

QpQr 


for  0  <  a:  <  Xj 


/(:(0=sw-x)= 


MiniQp,Qji) 

QpQr 


for  x^<x<X2 


5Ms£. _ £ 

QpQr 


for  X2<X<XMax 


0  otherwise; 


(4.2) 


and 


Xi  =MiniQp,Qg); 

^Max  ~  Qp  Qr'  (4.3) 

^2  “  ^Max  ”■  • 


It  is  easy  to  see  that  /(.)  >  0  for  all  x .  It  is  also  true  that 


^Max 

Jfa(t)  =  SW-x)dx  =  l. 


To  show  this,  we  begin  by  integrating  the  first  term  of  f(.)  with  respect  to  x . 
”  1 


o(3f>QK  QpQr  2 


-dx=- 


x^  (Min(Qp,Qp))^ 


^QpQr 


^QpQr 


(4.4) 


Integrating  the  second  term  with  respect  to  x  gives 


Mm(QpM^._Min(Qp,Qp)  _  Mm(Qp,Q^)  _ 
I  QpQr  QpQr  QpQr  ^ 


MinjQp^Qp) 

QpQr 

Mm(Qp,QR) 

QpQr 


(^Max  2X|) 


[Qp  +  Qr-  ^MQp/Qr  )]  • 


(4.5) 
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Integrating  the  third  term  gives 

9  I 

f  XMax-X^^_XMa^  FAto 

X  QpQr  QpQr  ^2  ^QpQr  P2 

^Mgx^2  ,  ^2  _  ^ar  ^Ator^l  |  ^  (4.6) 

QpQr  QpQr  ^QfQr  "^QfQr  ^QfQr  QpQr  ^QpQr 
^  (xm^ -x?f  _  X?  _  (Min(0p,0v )f 
2QpQi?  2QpQr  “^QpQr 


Stimming  the  results  of  equations  (4.4),  (4.5),  and  (4.6)  gives 

~  +  Qr  ~  A4m(Qp, Qp )] . 

If  Min{Qp,Qp)  =  Qp  then 


If  Min{Qp,Qp)  =  Qp  then 


Therefore,  it  is  true  that 

^Ato 

J  /(f(0=SlV-x)da:=l. 

0 

As  a  consequence,  /(1(f)  =  SW  -  x)  satisfies  all  of  the  requirement  to  be  a  density 
function. 
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C.  The  Net  Inventory  Distribution 

Baker’s  equation  (38)  (page  30,  Reference.[l])  provides  the  structure  for 
this  approximation  to  that  Poisson  case.  We  will  assume  that  the  demand  during 
the  aggregate  lead  time  is  Normally  distributed  with  a  mean  of  /z  and  a  standard 
deviation  of  cr.  Now  if  this  is  to  be  an  approximation  to  the  Poisson  case  then 

<T  =  -JJi  should  be  true.  However,  there  are  items  managed  by  the  Navy's 
Inventory  Control  Point  which  don't  meet  that  condition.  Thus,  we  will  assume 
that  (T  can  take  on  any  value  but  that  will  be  sufficiently  large  that  it  will 

exceed  3ct  so  that  there  will  be  a  negligible  probability  of  a  negative  demand 
during  the  aggregate  lead  time. 

Finally,  we  will  follow  the  approach  of  Hadley  and  Whitin  [3]  in  assuming 
that  demand  is  truly  continuous  (see  Chapter  4  page  192)  rather  than  attempting 
to  deal  with  a  continuity  correction. 

We  first  rewrite  equation  (2.1)  which  is  Baker's  equation  (38)  as  follows: 
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forO<z<Xi 


k:  J 


QpQr  , 


;=0 


/[N(t)  =  SlV-z]  =  -^ 


QpQr 


jjf{z-j)dj 


QpQr 

+ 


1  ^ 


for  Xi<z<X2 


(4.7) 


for  z>X2 


We  begin  with  the  first  part  of  this  equation,  the  case  where  0  <  2  ^  rcj . 
Note  that,  as  in  Chapter  2,  z  r^resents  the  amount  that  tiie  net  inventory  is 
below  SW  at  any  given  time  t.  In  addition,  tmder  the  current  assumption  that 
demand  during  the  aggregate  lead  time  is  Normal,  it  follows  that  if  2  is  the 
demand  ttien 


/(x)=- 


1 


-i(^y 

^  a  J  . 


a-JlK 

where  n  is  the  mean  demand  during  the  aggregate  lead  time  and  cr  is  the 
standard  deviation.  The  first  line  of  equation  (4.7)  when  equation  (4.8)  is 
introduced  is 


(4.8) 


1 

QpQr 


1 

a-j2^QjQ 


(4.9) 
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The  right-hand  side  of  equation  (4.9)  tends  to  be  awkward  to  use,  especially  in 
developing  the  integral  upper  and  lower  bounds  when  converting  (4.8)  to  a 
standard  Normal  (i.e.,  N(0,1))  density  function.  Therefore,  let  v  =  z-j.  Then 
j  =  z-v  and  dj  =  -dv.  When  j=0,v  =  z  and  when  j  =  z,v  =0.  The  left-hand 
side  of  equation  (4.9)  changes  therefore  to 

Next  define 


r  =  - — —  and  0(r)=-7=e  2 


Then, 


f{v)  = 


(4.11) 


From  the  definition  of  r,  it  follows  that  dr  =  —  and  therefore  f  (v)dv  =  <|)(r)dr 

G 

since  the  o  terms  cancel. 

Next,  we  need  the  upper  and  lower  boimds  for  r  in  the  integral.  They  are 
obtained  from  substituting  the  bounds  for  v  into  the  equation  defining  r  above. 
When  r;  =  0  and  t?  =  z. 


r  =  ~  and  r  =  - — respectively 

<y  a 

Finally,  we  need  to  transform  z-v  into  a  fimction  of  r.  From  the  definition  of  r  it 
follows  that 

v  =  ^  +  rc  and  therefore  z-v  =  z- jn-ro . 

Putting  all  of  tiiis  information  together. 


58 


a 


=  J  {z-ii-rc)(l>{r)dr. 


QpQr 


The  first  part  of  the  right-hand  side  of  equation  (4.12)  is: 


(4.12) 


a 


G 


— ^  J  (2-ju)^(r)dr  =  -^-^  J  <lKr)dr  =  ^^ 

QpQr  ^  QpQr  ^  QpQr 


\  G 


\  O 


(4.13) 


where  0(m)  =  J  q>{r)dr .  The  second  part  of  the  right-hand  side  makes  use  of 


Identity  1  of  Appendix  4  of  Reference  [3]. 
2-u 


-1 


a 


z-M 

.  c  . 


(4.14) 


QpQr  ^  Q>Qr  ^  Q>Qr 

o  a 

Adding  equations  (4.13)  and  (4.14)  gives  the  final  form  of  equation  (4.12). 


z-n 

a 


^  j{z-v)f{v)dv  =  -^^  J  (z-  fi-r(7)(t>ir)dr 


QpQr 


JL 

G 


(4.15) 


_  z-ili 

fz-z^Y 

a 

/h\ 

fz-ju'l 

QpQr 

rUJ  ^ 

J 

QpQr 

<P\ 

1  <T  j_ 

For  the  case  when  Xj  <  z  <X2,  the  density  function  for  the  net  inventory  is 
(see  equation(4.7)) 


Following  the  same  arguments  as  for  the  derivation  leading  to  (4.15),  the 
final  form  for  the  first  part  is 
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(4.17) 


4> 


f  z  —  'r^-  —  11^  (  r  —  1/  ^ 


V  O 
a 


iZif 

V  cT  ; 


QpQr 


The  final  form  for  the  second  part  is 


QpQr 


QpQr 


\Jiij 


-a> 


JJi 


(4.18) 


Adding  equations  (4.17)  and  (4.18)  give  the  final  form  for  equation  (4.16). 


z-x^- 


2-M 

a 


QpQr 

_ ^ 


QpQr 


n-i 


Z-Xi- 


(4.19) 


For  the  case  where  X2  £  2,  the  last  part  of  equation  (4.7)  applies.  It  is; 


QpQr  5 


QpQr 


QpQr. 

MiniZrXu^) 

J  (XMox-Mz-M- 


(4.20) 


In  this  case,  the  first  term  is  a  repeat  of  the  first  part  of  equation  (4.16).  The 
second  term  is  similar  to  the  last  part  of  equation  (4.19)  but  its  upper  bound 
changes  now  to  X2-  Therefore,  their  sum  in  final  form  is 


QpQr 


QpQr 


o 


-0 


£-£ 

.  (T  . 


(4.21) 


QpQr 
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The  remaining  term  of  (4.20),  upon  conversion  of  j  to  z-v,  is 


J  /(^)^^+7^7r-  J  vfiv)dv, 

z-Min(z,x^)  ^P^Rz-Mm,x^) 


which,  when  v  =  n  +  ra,  becomes 


^0,0  J  r<(>(r)dr. 

^^R  z~Min(z,XM^)-n  z-Min(z,XM^)-n 


This  then  reduces  to  the  following  final  form: 


QpQr 


+  „g 

QpQr 


z~Miriz,x^ 


(4.22) 


Summing  equations  (4.21)  and  (4.22)  provides  the  final  form  for  equation 


(4.20). 


QpQr 


J  ;/(2  -j)dj  +  J f(z-j)dj 


Min(z,x^) 


=  ^2LiL  $ 

QpQr  L 


QpQr  L  v  ^  /  v  cr 


■K^ 


(4.23) 


QpQrL  ct  j  ct  ;j 
^x^-z+n\Jz-Min(z,XMa:c)-B\  Jz-X2-^l 
QpQr  L  V  o  J  \  (7 

+  _2_r/ Z-Mm(z,XMax)-B]  /  2-^2-MY 
QpQrLV  <T  j  C7  Jj' 
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The  final  components  of  the  density  function  for  the  net  inventory  are 
equations  (4.15),  (4.19),  and  (4.23).  However,  to  reduce  the  visual  complexity  of 
the  density  function  for  the  net  inventory  we  will  use  a  shorthand  notation; 
namely. 


-  x-^ 

X  = — — 


a 


(4.24) 


With  this  notation  the  density  function  for  the  net  inventory  assuming  a  Normal 
distribution  for  the  demand  during  the  aggregate  lead  time  can  be  written  as: 


QpQr 


[o(5)-<I>(2)]-^~~[^6)-0(2)]  for  0<z  <Xi; 


QpQr 


/[N(f)=SW-2]  = 


QpQr 


[^>(2 -Xj  ) -$(2)]  -  ^^-[0(2 -Xj )  -  0(2)] 

[4>(5)-0(2-Xi)] 


for  Xi<Z  <22,* 


^[«2-x-i)-<lx'z)]-^[4,(z-xO-0(h] 


QpQr 

, 


QpQr 


[a>(2-X2)-^»(2-Xi)] 


QpQr 

a 


+£M^__£±£J,j>C2  _  Miii(z/x^))-^(z  -X2)] 

[0(2 - Min( z,xj^ )) - 0( 2-X2 )]  for  X2<z. 

(4.25) 


.  QpQr 


Finally,  as  we  did  with  equation  (2.3),  we  will  rewrite  equation  (4.25)  with 
ths  QpQr  term  removed  to  simplify  the  derivation  of  Pquj{SW). 
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(z- /i)[o(6 j -<E>(2)] - cy[^(0) - (f>(zj^  for  0<z< Xy 


(z  -  z  -Xi )  -<E>fz)]  -  a[<l>(  z  -Xi)  -  z)] 

+Xi[4)(0)-^('z-Xi)]  for  Xi<z<x2; 


QpQRf[N(t)=SW-z]=\ 


(z-fi)\fo(z  -Xi)-^(z)]-(T[<t>(z-Xi)- <(f(h] 
+Xi[^(z^X2)  -0(2^;] 

+(XMax  -Z  +  n{^(z-Miliz^Max)) -0( 7^)] 
+<j[(l>(z -  Min(z,Xf^x))-(l>(z-X2)]  for  X2<z. 


(4.26) 


D.  Derivation  of  PoujfSW) 

1.  The  case  of  0  ^  SW  ^  jcj 

We  begin  by  writing  PQUjiSW)  in  the  structure  of  equation  (2.12). 

fi)iir(^^)  =  J  {(z-iu)[<J>(9)-0(2)]-cy[^5)-^(z)]|^iz 

vJpWr  ^  ^ 

(z  -  /f)[<E>(z-Xi)  -  <I>(2)]  -  o[^{z-Xi)  -  0{z)]j 

+xi[a)(6) -$(“)] 

(z - ^)[^(2~Xi)- ^(z)] -  o[^(z^^) - 0(2)] 
+Zl[^>(2^)  -  a>(2^)] 

+(^Max  - Z  +  M)[^(z-Mm(2,Zjv|^)) - 0(2^2)] 
+(T[0(2-Mm(z,XAto))  -  0(2 -3:2)] 


QpQr 


"•Z 

J 


dz  (4.27) 


QpQr 


J 


kz. 


As  we  did  in  Chapter  2  to  get  equation  (2.13),  we  collect  common  terms  and 
rewrite  equation  (4.27)  as: 
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QpQrPout(^^)='  J  (z- pL)^itl)dz-  J  (z- ^L)^(z)dz-a  j  0(5)dz 
sw  sw  sw 

CO  oo  oo 

+  (tJ  0(z)rfz  + J(z-/l)^>(?^)dz-crj0(2’^)rfz 

sw  X,  X, 

X2  CO  00 

+  Xj  J  <I>(0)rfz  -  J  0(z  -Xi)dz + I  ^{z-X2)dz 

X,  a:j 

^Mmx  ^Aiu 

+  /  {XMax-Z+ ft)^(^)dz+a  j  0ib)dz 

X2  X2 

CO  00 

-  J  (^Max  -  2  +  fl)Mz^)dz  -  (T I  (f>{z-X2)dz 

00  00 


(4.28) 


We  will  consider  the  terms  in  equation  (4.28)  which  contain  either  <l>(5) 
and  0(6)  as  constants.  We  will  evaluate  them  first. 


j  (z-/i)4>(5)dz=4>(6)]-y-/iz 
sw 


w 


=  0(5)  ^ 


■SW^ 


-/i(xi-SW)  ; 


-a  J  0(6)rf2  =  -o-0(6)[xi  -  SW] ; 
sw 


^2  ^ 

Xi  J  <I>(5)dz  =  2iO(6)[a:2  -  2i]  =  0(6)[2iX2  -xf] , 

In  the  next  term  we  make  use  of  the  definition  of  X2  given  by  equation 
(4.3).  We  rewrite  it  as 


■"  ^Max  ^2* 
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^Max 

J  (Xmox  -Z  +  ^lM^)dz  =  {Xmox  +  “ ^2] " 


'1  y.2 

^Max  ^2 


=  +  ii)x^  -  +^?) 


=  0(5) 
=  0(6) 
=  0(5) 


/ai  + 


£M3£i_££^ 


=0(5) 


iU^l  + 


(^M«x-^2)^l 


a  J  0(5)d2  =  cT0(5)[xjvto  - ^2] = <^<K^)xi . 


Smnrning  all  of  these  "constants"  gives 


-  /2(xi  -SW)\-  (T0(5)[xi  -  sw] 


I  .  2 

0(6)[a:i3C2-^i]  + 0(5)1 


Aai  +  - 


:c? 


+cr0(5)xi 


=  0(5)1  ———  +  /iSW  +  a!:iX2|  +  <70(5)5  W . 


Next,  recalling  equation(4.24),  we  know  that 

z-ti 


and,  therefore. 


z  =  - 


J*  dz  J  -pTr-r  SW-fi 

dz  = — ,  and  SW  = - — 

o  <7 


(4.29) 


As  a  consequence,  we  write  tiie  second  term  of  ecjuation  (4.28)  as  follows: 


CO  ®®  /  \ 

-  f  (2  -  ju)0(2)d2  =  -<7  f  I  ^  lo(2)d2  =  -o^  f  ^(2)<i2  =  -a V(^) ;  (4.30) 
sw  sw^  ^  ^ 
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where,  using  Identity  7  from  Appendix  4  of  Reference  [3],  we  define 

oo 

/(m)  =  J  r<E>(r)df  =  - 1(1-  m^)4>(m)  +  u0(m)|  . 

U 

The  fourth  term  of  equation  (4.28)  can  be  easily  reduced. 

oo  oo 

(tJ  ^z)dz  =  a^  j  0{z)d2  =  tT^<I>(^).  (4.31) 

sw 

oo 

To  reduce  J (z  -  fi)^{z-Xi)dz  we  need  to  introduce  a  change  of  variable. 

Let  M  =  z—  Xj .  Then  u  - - ^ — — ,  and  therefore  z — /£  =  Xj  +  gu  .  Finally,  when 

z  =  Zj ,  u  =  6 .  We  also  need  to  add  a  definition  for  odu).  It  is 

oo 

a{u)  =  J  <I>(r  )df  =  (piu)  (4.32) 

U 

Therefore,  we  can  write 

oo  oo 

I  (z  -ju)^>(z-Xi)dz  =  crj  (xj  +Gu)^{u)du 

a:,  6 

oo  oo 

=  cXTi  J  <E>(m  )du+cf^j  u^(u)du  (4.33) 

6  t) 

=  axja(b)+G^Y(b). 

oo 

Next,  we  need  to  reduce  -gJ  (j){z-Xi)dz.  We  get 

oo  oo 

-CT  J  ^(z-x{)dz  =  -CF^  j  =  -ct^4)(6)  .  (4.34) 

a:,  D 

The  eighth  and  ninth  terms  from  equation  (4.28)  cancel  as  shown  below.  This  is  a 
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consequence  of  letting  u  =  ^  ^  in  the  eighth  term  and  then  u  =  - — ^ — — 

cr  a 

in  the  ninth  term. 

oo  oo  oo  oo 

-Xi  J 0(z-  X|)dz  J  -  X2)dz  =  -Xioj ^iu)du+Xi<jj  ^^)du=  0 . 

The  last  two  lines  of  terms  in  equation  (4.28)  simplify  down  to  one  term 
because  the  rest  of  the  reduced  terms  cancel.  As  in  the  preceding  derivation,  u  is 
defined  differently  for  each  part.  First, 


-J  (Xmox -z  +  ij)^z-X2)dz  =  -<7Xjf^{^(u)du+ax2 

^2  U 


^u)du  +  <J^  [ u^(u)du 


6 


=  +  OX20it))  +  <7^(5) 

=  -OXia{h)  +  <T  V(5) . 


OO  eo  oo  oo 

J  (^Ate  - 2  +  XMax)dz  =  OXMaxj ^(u)du  -  OX^^xj  <P(M)du-G^j tMu)du 

6  5  6 

=  oxMaxCdb)  -  atAtoa(5)  -  aV(9) 

=  -cr^r(O) . 

FinaUy, 


oo  oo  oo  oo 

-CF  J0(z -X2)dz+a  j  <t>{z-Xj^^)dz  =  -G^  j<l>{u)du+a^j(fi{u)du=  0. 

xj  Xito  0  6 


When  the  terms  in  the  last  two  lines  of  equation  (4.28)  are  summed  the  result  is 
merely  -<xcia(5). 


67 


-(^Max-Z  +  I^)j  ^(Z-X2)dz-CT  j<l>(z-X2)dz 

X,  Xj 

oo  oo 

+  (XMax-Z+  j  ^Z-XMax)d^  +  <^  J  ‘t>(z‘-xi^)dz  (4.35) 

^Max 

=  -oxia(b)  +  G^Y(t)) -o^y(^)  =  -mjafO) . 

Sumnung  the  results  of  equations(4.29),  (4.30),  (4.31),  (4.33),  (4.34),  and 
(4.35)  gives  the  final  form  for  Pout(^'^)  of  0  <  SW  <  x^. 

Pout  (SW)  =  0(5)|-^ + iiSW  +  XiX2  > + oSW^(5) 

-  (T^x(SW) + <T^C>(SW) + <xtia(6) 

+  crV(&)  -  <t^O(D)  -  <jria(5)  +  (tV(5)  -  cr^y(5)  (4.36) 

=  4>(5)|-^^  +  nSW+  X1X2  j  +  aSWpib) 

-  (7^[y(^)  -  ^(SV^j + <T^[y(5)  -  0(5)] . 

2.  The  case  of  Xi  <  SW  <  JC2 

In  this  case  the  first  part  of  equation  (4.27)  disappears  and  the  second  part 
has  SW  as  the  lower  bound  on  the  integrals.  The  initial  form  of  Pout(^^)  is 
now 

I  (z-^)[^(z-Xi)-^h]-o[<l>(z-Xi)-(l>(z)] 
ibirr(SW)  =  ^  ■  1  ^  ——1 

QpQr  5w  |^+Xi[O(0)  -0(2  -Xj  )] 


,W2.(4.37) 

>] 

> 


(z  -  jtx)[o(z  -x^ )  -  0(z)]  -  c[<l)(z  -Xi )  -  0(2)] 

I  “  +Zi[0(z  -Z2)-0(2  ^l)] 

. . — '  I 

QpQr  +(XMax  -2  +  -  Miriz^^ax'»  -^(z  -^2 

+cr[0(z  -  Mfn(2,Zj^^ ))  -  0(2 -Z2)] 
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Collecting  like  terms,  we  write  the  equivalent  form  to  equation  (4.28)  for 
this  case. 


OO  oo 

QpQrPoijt(SW)=-  j  (z-ii)4>(z)dz+G  J  <t>Cz)dz 
sw  sw 

oo  oo 

+  J  (z  -  /i)<I>(z  -Xi)dz  -  a  J  ^(z^Xi)dz 

sw  sw 

oo  oo 

+  Zi  J  0(5)iiz-zi  J  4>(z  -  Xi)dz  +  Xj  J  <^(z-X2)dz 

sw  sw 

^Mex  ^Max 

+  (Zjvtax  ~  2  +  /^)  J  +<J  j  ^(5)dz 

oo  oo 

~  -2  +/i)  J ^(z-X2)dz-a  j  <t>{z-X2)dz 

oo  oo 

~z +/i)  J  4>(z — J  <l){z  — x^^fj^)dz. 


(4.38) 


The  last  three  lines  of  equation  (4.38)  are  identical  to  the  last  three  lines  of 
equation  (4.28).  In  addition,  the  first  line  of  equation  (4.38)  is  the  same  as  the 
second  and  fourth  terms  of  equation  (4.28).  Therefore,  we  need  only  derive  the 
reduced  forms  for  the  second  and  tiiird  lines  of  equation  (4.38). 

For  the  terms  of  the  second  line,  let  m  =  z  -  Xj .  Then,  when 

z  =  SW,  u  =  SW  -  Xi .  The  rest  follows  from  the  definition  of  u  above  (just 
preceding  the  definition  of  a(u)  given  by  equation  (4.32)). 

oo  oo 

J  (z  - ju)0(z  -Xi)dz  =  a  J  ixi  +  ou)^{u)du 
sw  sw-:r,  (4.39) 

=  CTXia(SW-Zi)  +(T^y(SW-Zi) . 
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oo 


oo 

-(tJ  ^(z  ^Xi)dz  =  -(p-  J  0(u)rfn  =  -g^O(S W  -Xi),  (4.40) 

sw  snfe; 

The  first  term  of  the  third  line  is  a  constant  term; 

Xi  J  <I>(5)d2  =  ri0(5)(r2-SW).  (4.41) 

sw 

The  second  term  of  the  fliird  line  is 


-Xi  J  0(2  -  Xi)dz  =  -oxi  J  ^{u)du  =  -(XCia(SW  -  Zi) .  (4.42) 

sw 

The  third  term  is  the  same  as  for  the  previous  case.  There  it  was  not  completely 
reduced  since  it  cancelled  the  second  term.  Here,  it  needs  to  be  reduced  since  its 
lower  botmd  is  different  from  the  second  term. 


00  00 

Xj  J  0(2  -  X2)dz  =  Xioj  ^{u)du  =  Xiacx(b) .  (4.43) 

X2  0 


Therefore,  adding  die  results  of  equations  (4.39)  through  (4.43)  gives 

aria(SW-Xi  )+ a^y{SW  -x])  -  g^^{SW  -  Xj) 

+  Xi^>(0)(x2  -  SW)  -  oxia(SW  -Xj) + XiOor(S) 


=  <T  V(5W  -  xj)  -  ct^4>(SW  -  xj) 
+  Xi^(5)(x2  -  SW)  +  Xi<7a(5) 


(4.44) 


Combining  this  result  with  the  unchanged  derivations  for  the  previous 
case,  we  get 

+  ju-Sw|+ axj0(6) 

-  (T^[r(W)  -  ^{m)]  (4.45) 

+  <T^[y(SV\^Xi)  -  <1>(SW-Xi)] . 
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3.  The  case  of  X2  ^  SW  <  x^ax 


Equation  (4.37)  can  be  easily  modified  to  get  the  starting  form  of 
Pou7’(SW)  for  this  case. 


Pout(^^)  = 


QpQr 


I 

SW 


(z-  n^^iz-Xi)  -  0(z)]  -  ay>iz  ^Xi)  -  0(z)] 
+Zi[4>(r^)  -  O  (z^)] 

•^Xmox  - 2  +  tl\^(.Z-Miniz,XMax)) -^(z^X2 )] 

+a[(l>{z-  Min{z,xj^))  -  ^(z^2)] 


^.(4.46) 


Because  of  the  Mm(z,z^^)  term  it  is  appropriate  to  consider  the  interval 
between  X2  and  by  itself.  Equation  (4.46)  can  be  rewritten  as 


00  00 

QpQrPout(^^)  =  -  J  {z-n)<M,z)dz+a  J  0(z)dz 

SW  SW 

00  00 

+  J  (z  -  //)0(z  -Xi)dz  - o  J  (l>(z~  Xi)dz 

SW  SW 

00  CO 

J  ^z-Xi)dz+Xi  J  ^(z-X2)dz 


SW 


SW 


(^Moc-z+Ai)  J  <I>(5)dz+cT  J  <l>{b)dz 


SW 


SW 


(4.47) 


00  CO 

-(xmox-z  +  h)  f  <^{z^X2)dz-a  f  0(z -X2)dz 

SW  SW 

00  CO 

+  (2Mfl;c-2  +  Ai)  J  <^{z-X}^)dz  +  o  J  <j>{z-~Xj^^^)dz. 


'^Max 


The  first,  second,  and  last  lines  are  the  same  as  for  the  previous  case.  The 
first  term  of  the  third  line  does  not  change  either.  The  reduction  of  the  remaining 
terms  follows. 
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Xi  J  <E>(z  -  X2  )dz  =  aX|  J  <^{u)du  =  oxia(SW -  X2 ) .  (4.48) 

SW 


J  (^Max  -  2+  /X)<K5M2  =  iX}^^  +  /i)<I>(5)(2^ax  -5W)  -4)(6)-~b 

SW  ^  P 


Max 

isw 


A  A  r*- 

=  r^M<rc  +  MmOXXMax  SW)-  <^0)  ^  -  ^ 


(4.49) 


=  4>(6) 


j(XMaxSW)^+fl(XMax-SW) 


^Max 

a  J  <l>{t))dz  =  cr0(5)(xAto  -  SW) .  (4.50) 

SW 


-  I  (2A4ax-2  +  /i)^(2-X2)fi2 
SW 

oo  oo  oo 

=  -^Max<^  J  ^(u)du  +  0X2  J  ^>(M)d«  +  <T^  J  U^{u)du 

SW^x^  sw-x^  sW-x^ 

=  -^Mflx<«3J(SW^X2)  +oa:2a(SW"-X2)  +  <tV(SW^)  (4.51) 


=  -Xi(Ta(SW  -  X2)  +  a^yiSW  -X2). 


-a  J  0(2  -  X'^dz  =  -o^  J  ^(Si)du  =  -<t^O(SIV  -  X2 ) .  (4.52) 

SW  5w-Xj 

Summing  up  the  resialts  of  equations  (4.48)  through(4.52), 

aX^(SV^-X2)  +a>(t))|i^  +^l(XMaxSW)^ 

+  (T(l>(t))(xi^„x  -  5W)  -  Xi<ja(SW-  X2)  +  a^(SW  -  X2)  -  ct^4>(SW  -  X2)  ^  ^  ^ 
=  0(t))|--^‘^^  +)l(XMax  SW)^  +  G$(t))(XMax -SW) 

+  a'\y(SW- X2)  -  a>(SW-X2)] . 
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Summing  the  other  terms  which  were  unchanged  from  the  previous  case, 

we  get 

axia{SW-Xi ) + cr^[y(SW^^)  -  4>(SW-Xi )] 

-  cT^[y(^)  -  a>(^)]  -  Gx^oiSW-x{)  -  <t^[7(5)  -  a>(5)] 

=  <T^[y(SW^)-<I>(SW^)]-  -4>(^] 

-(T"[y(6)-4»(5)]. 

Finally,  adding  equations  (4.53)  and  (4.54)  gives  final  form  for  PoUT^^^ 
for  the  case  of  X2  <  SW  <  • 

PoiniSW)  =  4>(5).  +/((^Ator  -  SW)  ■  +  -  SW) 

-  CT^[y(^)  -  4»(^)] +(T^[r(SW-Xi )  -  <I>(SW-;ci)]  (4.55) 

+  cr2[y(5w^)-^(5w^)]-  (T2[7(a)-<D(5)]. 

4.  The  case  of  SW  >  Xj^^x- 

This  case  is  similar  to  the  preceding  one  in  the  fact  that  equation  (4.46) 
also  applies  here.  Now,  however,  Min{z,Xn^)  =  Xj^.  Thus,  equation  (4.46) 
takes  on  the  following  form: 

(2  -  ju)[4>(2~2i)  -  ^(2)]  -  <y[0(2-2i)  -  0(2 
+Xi[<E>(2  -22)  -  <I>(2  -2i)j 
"^^Max  “  2  A^)^*J*(2  —  ^Max^  ~  ^(-2  ~  ^2)] 

+(^(l>{z-  XMax)  -  0(2^22)] 

The  fourth  line  of  equation  (4.47)  is  now  eliminated.  The  reduced  version  of 
equation  (4.47)  is  given  by  equation  (4.57)  below.  Only  the  terms  in  the  last  line 
of  equation  (4.57)  have  not  been  analyzed  somewhere  above. 
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oo  OO 

QpQrPout(^^)  =  -  J  {z-fl)^{z)dz+C7  J  (}>Cz)dz 
sw  sw 

OO  OO 

+  J  (2  -  jU)<I>(2  -  Xi)dz  -  or  J  ^{z-  Xi)dz 
sw  sw 

OO  OO 

J  ^iz-Xi)dz  +  Xi  J  ^>(2- 22)^2  (4.57) 

sw  sw 

OO  OO 

-  J  (XMax-Z+fl)^Z^2)dz-G  j  <t>iz-X2)dz 


SW 


sw 


J  ^Max  ^  ^)^iz  Xf^^)dz  j  0(2  Z^^)dz 

SW  sw 

The  tenns  of  the  last  line  of  equation  (4.57)  are 


J  iXMax-Z+fiMz-XMax)dz=-a  j  ^\^{z  - XMax)<iz 

sw  sw ^ 


=  -cy^  ]  u^{u)du  =  -crV(SW  -  x^ax ) , 


(4.58) 


SW-X, 


and 


o  j  <li{z-XMax)dz=^cy^  J  (l)Cu)d7i=c^^iSW-XMax)-  (4-59) 


SW 


The  sum  of  equations  (4.58)  and  (4.59)  is 


J  (xmox  -  2+  /i)0(2-2Aia^)d2+cr  J  <liz-Xj^)dz 


sw 


sw 


(4.60) 


=  -(T^[r(sw  -xm„)-  «(sw  -  im..)]  ■ 


The  final  form  for  Pojjj(SW)  for  tihe  case  oixj^<  SW  is,  therefore. 
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(4.61) 


Pout(sw) = -<T2[r(SW)  -  4KSW)] 

+  (T^[y(SW  -Xj)-^SlV  -A^i)] 

+  <72[y(SW^)-®(SW^)] 

-  (T^IrCSW  -  Xjto)  -  4KSW  -x„„)] . 


E.  Summary 

This  completes  ttie  derivations  of  the  equations  for  Pqijj’{SW)  imder  the 
assumption  that  demand  during  the  aggregate  lead  time  is  Normally  distributed. 
For  0  <  SW  <  Xi,  equation  (4.36)  applies;  for  Xi  <  SW  <  ^2,  equation  (4.45)  applies; 
for  X2  S  SW  <  Xj^,  equation  (4.55)  applies;  and,  for  Xj^g^  <SW,  equation  (4.61) 
applies. 
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CHAPTER  5  -  DERIVATION  OF  B(SW)  FOR  THE  NORMAL 

MODIFICATION 


A.  Introduction 

This  chapter  presents  the  derivations  of  the  formxilas  for  all  of  the  cases 
for  the  expected  number  of  backorders  at  any  instant  of  time,  B{SW),  for  the 
Baker  model  where  demand  during  aggregate  lead  time  is  assumed  to  be 
Normally  distributed.  The  density  function  for  the  Normal  modification  to  ttie 
Baker  model  was  developed  in  Chapter  4  and  is  given  by  equation  (4.25).  The 
derivations  in  this  chapter  will  also  rely  heavily  on  the  derivations  in  Chapter  3 
which  gave  the  formulas  for  B(SW)  for  the  Baker  model  [1).  In  place  of  equation 
(3.1)  we  now  have 


oo 

B(SW)=  J  {z-SW)J[m)=SW-zYz  (5.1) 

sw 

Next,  we  can  separate  this  equation  into  the  following  two  parts; 

oo  oo 

B(SW)=  J  zJ[N{t)  =  SW-z]dz-SW  j  J[N{t)=  SW- z]iz . 
sw  sw 

As  we  noted  in  Chapter  3,  the  second  term  is  nothing  more  tiian  the  product 
SW PQUjiSW).  The  derivations  for  Pqijj(SW )  were  presented  in  Chapter  4  for 
the  Normal  modification.  Therefore  we  can  focus  our  attention  on  the 
derivations  associated  with  the  first  part.  We  will  denote  this  part  as 

oo 

Bi(SW)  =  j  zf[N(t)  =  SW  -  2]d2 .  (5.2) 

sw 

Therefore,  we  can  write 

B(SW)  =  Bi(SW) -SWPomiSW) .  (5.3) 
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B.  Derivation  of  BiiSW) 


1.  The  case  of  0  <  S  W  <  Xi 

For  the  case  of  0  <  S  W  <  x^,  when  the  net  inventory  probability 
distribution  is  given  by  equation  (4.25)  we  can  write  Bi(SW)  as 

B,(SIV)  =  .^  I  z{(2-A)[®(6)-«(j)]-<T[^(6)-^(2)]}lz 


sw 


1  7 


(z  -  /i)[^(2  -Xi)-  $(z)]  - (y[0(2  -  2i)  -  ^(z)]] 

+2i[0(6)-  ^>(2^)] 


fdz 


—I 

.r)„  J 


QpQr 


z< 


QpQr 


[0(2  -  Zj)  -  0(2)]-  ^^[0(2~2i)  -  0(2)] 


+  -  Al'V2,x„„))-^(2“jr2)] 


k 


Next,  we  rewrite  equation  (5.4)  into  its  individual  component  integrals. 


(5.4) 
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oo 


QpQRBiiSW)=  I  z(z-/i)0(D)rf2-  J  z(z-;i)a>(z)iiz-(j0(0)  J  zdz 
sw  sw  sw 

oo  oo  OO 

+  <7  J  z<}>(z)dz  +  J z(z  -  ju)<>(z  -Xi)dz  -  <jj  z^z-Xi)dz 

X2  00  00 

ZiO(6)  J  zdz  -  Xj  J  z^{z-Xi)dz + Xj  J  z^(z^l^)dz 


+  . 


+  J  -  z  +  -  Min(z,  x^^  ))dz 

00 

-  J  Zix^ax  “2  +  n)<Hz-  X2)dz 

X2 

OO  00 

+  or  J  z<^z-Min{z,Xi^))dz  -  <t  J  z0(2^^)dz . 


(5.5) 


The  first  term  of  equation  (5.5)  is  a  constant 

^  XX 

J  z(z-ju)0(5)iz  =  <I>(5)  j  z^dz-/x^(5)  J  zdz 


sw 


sw 


sw 


=  0(5) 

=  0(5) 


SW^ 


-AiO(5M 


2  sw^ 


3  3  2  2 


2  2 
2' 


(5.6) 


For  the  second  term  of  equation  (5.5)  we  recall  from  Chapter  4  that 

2  =  ^  ^ ;  sothat  z=az+fi,  z- fi  =  az,  and  dz=^cdz. 
o 

Therefore, 
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(5.7) 


-  {  z(z-^i)^(z)dz  =  -<7'^  {\cjz^  +  HZ  j^(z)dz 

SW  sfc*-  ■* 

<» 

2^<I>('z)d2-<7^  J  24>(2)^iz, 

ST^ 

We  now  define  5(m)  by  using  Identity  8  from  Appendix  4  of  Reference  [3]. 


5(u)  =  j  /^(r)dr  +2  ^fw) - m «)| ,  (5.8) 


Thus, 


-  J  2('2-/l)0('2)dz  =  -<T^  f  Z^4>('2)^-(T^/i  f  24>(z)rf2 

SW  ^  SW 

= -(T^^C^)  -  o^/iyfSW) , 
where  y(.)  was  defined  in  Chapter  4. 

The  third  term  of  equation  (5.5)  is  a  constant. 


(5.9) 


-(700)  (  zdz=-O00) 

SW 


xj  SW^ 
2  2 


(5.10) 


The  fourth  term  is 

oo  oo 

G  f  Z0{z)dz=G^  f  iH+GZ)(l)(z)dz 
SW  SW  (5.11) 

=  /iCT^4>(5W)+<T^0(SW) . 

Here  we  made  use  of  Identity  1  of  Appendix  4  of  Reference  [3];  namely. 


Jf0(r)fir=  0(u). 

U 


To  reduce  the  fifth  terai  of  equation  (5.5), 
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jziz-iii)0{z-xi)dz , 

X: 

we  define 

*  Z-Xi-fl 
u=  * — — . 

a 

Then,  «  =  D  when  z  =  Xi,  and 

z  =  Xi+  pL  +  ou; 

Z- jU  =  Xj  +  (TM  ; 

so  that 

Z(z-fl)-(Xi  +CTM)('Xi+/1  +  (TM) 

=Xiixi-\-  pO  +a(2xi+  pL)u+  . 

Therefore, 

oo 

J z(z-  pL)^(z-  Xi)dz 

oo  oo 

=  GXi(xi  +  ii)\  u)du  +  <T^  (2xi  +  jti)  f  u^(  u)du  +a^\u%  (u)du  (5.12) 

ft  ft  ft 

=  axi(x^  +  ii)a(b)+o^  (2xi  +  iJ,)y(^)  +cr^5(6) . 

The  sixth  term  uses  the  same  definition  for  u  as  equation  (5.12)  did. 

oo  oo  oo 

-a{  z/^  -Xj)dz  =  (i,  +  ii)i<Hu)du-<j^t  uip(.u)du 

i,  h  4  (5.13) 

=  — <T^(Xj  +  /i)4>(0)  — <T^0(O) . 


The  seventh  term  is  a  constant  term. 
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(5.14) 


zdz  = 


2  2 


The  eighth  term  also  uses  the  same  identity  for  u  as  equation  (5.12)  did. 

oo  OO  oo 

-Xi  J  z<E>(2  -  Xi)dz  =  -a:i(x2  +  m)  J  Mu)dz-  axj  J  7^(u)dz 

oo  oo 

=  -oxi  (xj  +  ju)J  <S>(^)du  -  G^Xi  J  u^(u)du  (5.15) 

1)  0 

=  -oxi  {x-i  +  n)(x{Q)-  ) . 


For  the  ninth  term  we  define  u  as 


-  Z-Xo-U 
U-— . ^ — — 

a 


Then,  we  get 


Xi  J  Z<^{z-X2)dz  =  OXij  {X2+^l  +  ou)<>{u)du 
^2  i) 

oo  oo 

=  ax^ix2+  fi)j  4>(«)(iM +g^xA  iii<P(u)du  (5.16) 

6  T) 

=  aici(x2  +  fi)a(^)+GhiY(^) . 

Notice  that  when  we  add  equations  (5.15)  and  (5.16)  the  sum  is  zero  since  the 
terms  in  one  equation's  reduced  form  cancel  those  in  the  other. 


The  tenth  term  uses 
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G 

for  the  second  part  (the  integral  from  Xj^  to  oo).  The  derivation  of  the  reduced 
form  of  the  tenth  term  is: 

oo 

J  ^(^Max  -2  +  /i}4>(z -Min(z,xj^))dz 

=  J  z(XMa^-Z  +  H)^0)dz 

Xj 

OO 

+  I  ^(XMax  -  2  +  -  XMax)dz 

^Max 

^Max  ^Max 

=  -  J  z^^(0)dz  +(XMax+f^)j  Z^0)dz  (5.17) 


(XMox"^  +gu)^(u)£i 
D 


=  <E>(5)j 

-  ^^(xmox  +  (5)  -G^5{b). 


XMox  X2 
2  2 


The  eleventh  term  uses  the  same  definition  for  u  as  equation  (5.16). 

00 

-  J  xiXj^ax  -Z+^)^{Z-  X2)dz 

Xj 

00 

=  -CT J  (a:2  +  jU  +  Gu){xj^  -X2-  Gu)^iu)du  . 


"0 


Expanding  the  integrand  factors,  we  get 
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(X2  +/l  +  (JuXXj^ax  -^2  -  <^) 

=  [-CtV  -a(X2  +H)u  +  G(XMax  “^2^ 

=  |^-(7^M^  -<j(X2  +  fX)u  +  GXiU  +  (X2  +  fi)xi^- 

Therefore, 

oo 

-  J  ^(^Max  -  2+  ll)Mz^X2)dz 

OO- 

=  -a J -G(X2-\-  pi)u+GXiU+(X2-\- n)xi 
=  <7^Jm  ^(u)du-  G^(Xi-X2-fi)j^(u)i^-aXi(X2+IJ.)\ ^u)du 

R  ft  ^ 


0  6 
-3e/ft\  _2/ 


=  CT^^fO)  -  G^(Xi  -X2-  x)y(^)  -  OXl(X2  +X)oc(h)  . 


The  twelfth  term  uses  the  same  definition  of  u  as  equation  (5.17) 

oo 

ctJ  20(2  -  Min(z,  Xj^  )dz 

xj 

^max  oo  _ 

=  a  J  20(tl)d2+C7  J  z^z-X}^ax)dz 


=  (T0(u;  — 


FMflx  2 


2h 


J  +/^  +  Gu)<l>(u)du 


=  CT0(5) 


+  +  /i)$(C)  +  C7^0(6)  . 


The  final  tenn  of  equation  (5.5)  is 


(5.18) 


(5.19) 
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oo 


-<7  J  20(z  -  X2)dz  =  {X2+fl  +  cu)<tl(u)du 
X2  T) 

=  -<P‘{X2  +  . 


(5.20) 


The  final  form  of  equation  (5.5)  involves  summing  all  of  the  reduced 
forms  contained  in  equations  (5.6),  and  (5.9)  through  (5.20). 


QpQMSW)^<Hti) 


xf  SW^  /g?  nSW^ 
3  3  2  2 

2 


-  <x^S  (ST^)  -  a^^yd^) 


-om^-  2 


+  /f(T^4>(SW)  +  (T^0(SW) 


+  axi  (xj  +  ju)a(5)  +  (2xi  +  /t  )y(5) + c  ^5(5) 

- cr^(xi  +  /i)^>(6) - <T^0(t)) + Xi4>(5) 

-  ax^ixi  +  jit)a(5)  -  (T^xiy(5)  +  0x1  (X2  +  M)a(0) 


(5.21) 


+  <7^Xiy(5)+0(5) 


-^^+^HxMax+^^-) 


^Max 


+  J“)y(S)  -  <t^5(S) +<7^5(5)  -  cHxi  -  X2  -  /t)y(5) 

+  <^^iXMax  +  f^)^(^) 


-  axi{x2  +  M)a(S)  +cr0(6) 

+  <T^^(5)  -  C7^  (X2  +  /i)<E>(S)  -  C7^0(5)  . 


2  2 


Equation  (5.21)  can  be  substantially  reduced  by  collecting  all  of  the  like 
terms.  Many  then  cancel  each  other.  The  resulting  equation  for  BjfSW)  for  the 
case  where  0  <  SW  <  Xj  is  therefore: 
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Bi(SW)  = 


<E>(0) 


QpQr 
.  oi^>(5) 


QpQr 


QpQr 


^  ^  sw"  ti4,„  p;j  fvc}  f£W^ 

6663  2  22  2 

^Max  _  .^2.  ^ 

2  2  2  2 

-a[5(5W)  -  0(^)]  -  /£[/ (5W)  -  <I>(SW )] 

+a[5(0)  -0(5)]  +  /i[y  (5)  -  <^(0)] 


(5.22) 


2.  The  case  of  Xi  <  S  W  <  a:2 

Equation  (5.4)  changes  to  the  following  fonn  for  this  case  since  its  first  line 
no  longer  applies. 


Bi(SW)  = 


QpQr 


^2 

sw 


(z  -#i)[<K2  -  *1)  -  4>(2)]  -  o[^(2  -  Zj)  -  ^  (z)]] 

[+z,[4>(5)-4K?^)] 


(5.23) 


-I 

.Dr.  J 


QpQr 


z< 


+  ^  [<K2  -  MiMzAMax))  -  ^■(Z  ^2)] 

+  Q^[<Kz-^iniz,XMax))-<l>iz~X2)] 


\dz. 


Therefore,  equation  (5.5)  now  changes  to  the  following  form.  Note  that  the 
lower  boimds  on  the  integrals  in  equation  (5.5)  have  been  replaced  by  SIV . 
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oo  oo 

QpQk®i(5^)  =  -  J  z{z-pi)^(z)dz  +  a  J  z<l>iz)dz 
sw  sw 

oo  oo 

+  I  z{z  -  /i)0(z  -Xi)dz  -<j  j  z<l>iz~Xi)dz 
sw  sw 

X2  oo  00 

+  aci<E>(5)  J  zdz-Xi  J  zO(z -Xi)dz+Xij zO(z -X2)dz 

sw  sw 

00 

+  j  ^(^Max  -  Z  +  fl)^Z  -  Mm(z,  X^ax  ))^^2 

ATj 

OO 

-  J  z(^Atoc  -  2  +  fl)^(z-X2)dz 

X2 

00  00 

+  <7  J  z0(z  -Mm(z,x^(,x))dz  -  c  J  z^{z  -X'^dz . 

x^  x^ 


(5.24) 


The  reduced  forms  for  the  first  and  last  three  lines  remain  the  same  as  for 
the  preceding  case.  Thus,  we  need  only  derive  the  reduced  forms  for  tiie  second 

and  third  lines.  Three  terms  in  the  two  lines  involve  u  =  (z  -Xi)  =  - — — * — — . 

<7 

Therefore,  in  these  terms  we  will  use 

z  =  Xi+^  +  Gu; 

Z- H  =  Xi  +  GU) 

and,  when  z  =  SW,  u  =  (S W -x^).  The  derivations  of  the  reduced  forms  for  the 
three  terms  are  shown  below. 
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(5.25) 


OO  _  oo 

J  z{z- Ii)^(z-Xi)dz  =  <y  f  (x^  +  /i+  au)(xi  +cTu)^(u)du 
sw  (SW^,) 


=  (TXi(Xi+fl)  ±  ^(u)du  +  <j^  (2xi  +  fi)  j  u^(u)du 
(SW-xJ  (SVV^) 


+  a' 


A  2  ^  ^ 

u  ^(u)du 


-x^) 


=  oxi(xi  +  fi)a(SW -Xi)  +  o  ilx^  +  ^i)y(SW  - Xi ) +ct  5(SW 


oo  oo 

-(J  j  z(l>(z  ^Xi)dz  =  -<y^  f  Ctj  +  iJ.+Gu)$iu)du 

SW  snfef; 

oo  oo 

=  -<y^(xi  +/i)  J  <l){u)du-a^  J  ut]>(^)du 

=  -CtHxi  +  fl)^iSW-Xi)  -C7^0(SW-Xi)  . 


oo  oo 

-Xi  j  z^{z-Xi)dz  =  -axi  J  {xi+fi  +  (m)4>{u)du 

sw  SW^, 

=  -<xci(xi  +ti)a{SW  -Xi)-(j^Xiy{SW  -rcj) . 


The  remaining  term  is  a  constant. 

Xi4>(0)  J  2rfz=Xi<E»(S) 


sw 


SW^ 
2  2 


Adding  all  of  the  reduced  terms  for  equation  (5.24)  gives 


Xi). 


(5.26) 


(5.27) 


(5.28) 
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QpQRB^iSW)  =  -(7^d(^)  -  <TVr(^) 

+  /XO-^0(^)  +<T^0(^) 

+  (XCiCxj  +  fi)(xiSW  -x^)  +  cx^iTxi  +  /x)y(SW-Xi) 

+  cr^5(SW-;ci)  -  (x^  +  /i)<I>(SW^) 

-  a^<l){SW-Xi)  +  Xi^CO) 

-  oxiixi  +  fi)a{SW  -  jci)  -  o^x{f{SW  -  x{)  +oxiix2  +  ^i)(x{^ 

.A] 

2 

-  -  0^80) + <t^5(5)  -  oHxi  -X2-  ju)y(5) 

-  OJCi (^2  +  fl)cc0)  +  a<l>0)  (^Afax  +  ^^)^0) 

+  o^<f)0)  -  cr^  {x2  +  n)^{t))  -  cy^<l>0) . 


+  a\Y0)  +  ^0)\-^+^+(xMax+^^)  ^ 


(5.29) 


Agaiiv  after  collecting  like  terms,  a  substantially  reduction  in  equation 
(5.29)  occurs.  Thus,  for  the  case  of  x^  <  SW  <  X2  we  get,  as  the  final  form  for 
Bi(SW), 

Bj(S  w) = -2^5L  + iiiMa.  _  id  _ 

QpQr  6622  2 

,  (^<l>0)\xhax  4 

QpQrL  2  2_ 

-<t[B(SW)-0(SW)] 

+<t[5(SW-Xi)  -  0(SW-Xi)] 

+iu[y(SW-Xi)  -  <I>(SW-Xi)] 

-Xi[y(5)-O(0)] 
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3.  The  case  of  X2  <  SW  <  x^ax 

Only  the  third  part  of  equation  (5.4)  applies  for  this  case. 

.  1  7 


B,(SW)  =  -^  I  z| 

QpQr  ^  I  ^Max 

QpQr 


^  ^OrO  ^  ^  ~  ^Max))  -  ^ (Z  -  3:2  )] 


iz.(5.31) 


+  q^[0(2  -  Min{z,X}^^^))  -  (piz  -  X2 )] 

As  a  consequence,  equation  (5.5)  is  reduced  to 

00  00 

QpQr®i(5^)  =  -  J  2(2  -  /i)<^(z)dz  +  a  j  z<l>{z)dz 

SW  SW 

00  00 

+  J  2(2  -  //)^(2  -xi)dz  -  a  j  z<p{z  -Xi)dz 

SW  SW 

00  00 

-2i  J  Z^{z-Xi)dz+Xi  J  20(2-22)^2 

SW  SW 

00 

+  J  -  2  +/f)0(2  -  Min(z,XM^))dz 

SW 

00 

“  /  ^(^Max  -  2  +  /x)0(2  -X2)dz 
SW 

®o  00 

+  (tJ  z<p(z-Min{z,X}^ax))dz-(T  j  z<p(z-x])dz. 


The  first  five  terms  will  not  change  for  this  case.  We  will  first  focus  on  the 
terms  involving  (2^2)  •  those  terms  we  define  «  as  «  =  (^^1  = 


(5.32) 


so  that  2  =  ^2  -/i  -  (T«. 
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(5.33) 


Xi  J  z<E>(z  -X2)dz  =  axi  J  (^2  +  ju  +  (Tu)<E>(M)dM 

sw 

=  axi  (^2  +  ix)oe{SW  -^2)+ (SW  -  ^2 ) . 

The  next  term, 

00  00 

-  J  2(JCa1«x  -  Z  +  n)^(z^)dz  =  -G  I  {X2+tl  +  Ou){XMax  “ ^2  “  pi)^iu)du  , 

SW  SW^ 

has  a  long  integrand  when  the  factor  terms  are  multiplied  together.  Doing  that 
first, 

(X2+fl  +  Ou)(Xj^  -X2  -  au)  =  Xi  (X2  +  tl)  +  G(Xi-X2  +  fl)u  -  . 

Therefore, 

oo 

-  J  -Z+fl)^(z^X2)dz 

SW 

OO 

=  -axi(x2+n)  J  ^(u)du 

SV^J 

oo 

-^2 -ju)  J  u^(u)du 
+  J  ^f^(u)du 

SW-x^ 

=  -axi(x2  +  fj.)a(SW  -  X2)  (5.34) 

-  o^(xi  -  X2  -  ii)y(SW  -X2) 

+  g^8(SW~-X2). 


The  final  term  involving  (z  -Z2)  is 
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(5.35) 


CO  oo 

-O  J  20(2  -X-^dz  =  -CT^  J  (^2  +/i  +  (TM)0(u)rfM 
SW  SVV-ATj 


=  -<J‘ 


{x2  +  iU)^>(SW  -X2)  -a^(fi(SW  -X2) . 


The  remaining  two  terms  need  to  be  subdivided  in  the  intervals 
SW^2^a:^  and  z>Xf^.  In  tiie  latter  range,  m=-^— 

(T 


\W 

=  J  z(Xji^ajc -Z+fl)^(^)dz+  J  z(XMax-Z+ li)^(z-Xj^^^)dz 


SW 


=  a>(0)  J  ^(Xj^ax  +/^)  -  j (Xmox  +M  +  OU)l^(u)£l 

SW  5 


=  W)\(XMax+l^) 


Xmox  SW^  xlf^  ^  SW" 


3 

“  *2 


-<7^(xj^  +fi)ju^u)du-o^ju  <^(u)du 

-  ^(5)[ 1 

1  6  2  3  2  J 

-  c7^  +/i>)y(9)  -  a^5(0) . 


(5.36) 


a  J  z^{z  -  Min(z,Xj^x)dz  =  00(6)  J  2^2 

SW  SW 

..2  /’T«  ,2 


+  0^ 


{X}^^  + ^  +  auyi>(u)du 


=  (T0(6) 


2  2 


+  (^Max  +  /i)^(^)  +  ^0(S)  . 


(5.37) 


Adding  up  all  of  the  reduced  terms  for  equation  (5.32)  gives 
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QpQrBi{SW)  =  -G^sw)  -  G^nriSW) + fiG^^iSW) + ct^^(SW) 

+  oxi(xi  +  fi)a{SW-Xi)  +  gH2xi  +  ^i)yiSW  -Xi) 

+  cr^5(SW  -  Xi)  -G^(xi  +  ju)<E>(SW  -Xj) 

-  <T^0(SW  -  arj)  -  axiixi  +  fi)aiSW  -  x-^-a^x^yiSW  -arj) 
+  <xci(x2  +  jU)a(SW-X2)  +  o^Xiy{SW  -X2) 

- 0%i(ac2  +  fi)a(SW -X2)- g^x-^  -x^- n)yiSW -X2) 


+  cy^5(SW  -  ^2)  +  ju)^(SW  -X2)  -g^^(SW-X2) 

■  0(5)|— +^l)SW^  I  f^Maxl 


(5.38) 


6  2  3 

V. 

-  +  /<)y(6)  -  <7^5(5) +01^(6) 

+  O^Umwc  +  /i)4'(0)  +  <7^«(6)  ■ 


J 

SW^ 


Collecting  common  terms  results  in  the  final  version  of  Bi(SW)  for  the 
case  of  X2  <  SW  <  x^^. 


Bi(SW)= 


<I»(0) 


QpQr 


^Max  ^  SW^  ,  ^l4iax  iXMax  +  tl)SW^ 
6  3  2  2 

q'0(5) 


QpQr 


'2  ^ 

•^Max  _  ji2 

2  2 


QpQr 


+c7[B(SW  -  ;ci)  -  0(SW  -  xi)] 
+ju[7(SW^)-  ^>(SW-  Xi)] 
+a[5(SW-X2)  -0(5W^)] 
+n[y{SVr-X2)  -  ^>(SW-X2)] 


(5.39) 
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4.  The  case  of  SW  >  x^ax 

Equations  (5.31)  and  (5.32)  apply  here  also.  The  part  that  needs  no  longer 
to  be  considered  is  where  Min{z,Xj^a^)  =  z  which  was  the  case  when  SW  <  Xj^g^ . 
Thus,  equation  (5.31)  can  now  be  written  as 


Bi(SW)= 


QpQr 


h 

SW 


Wz.(5.40) 


Therefore,  equation  (5.32)  can  be  rewritten  as 


QpQi?®i(SW)  =  - J  z(z-/x)0(z)dz  +  (T  J  z0(z)dz 

SW  SW 

oo  oo 

+  J  z(z  -  fi)4>(z  -Xi)dz  -  (T  J  z0(z  -Xi)dz 

SW  SW 


-Xi  j  zO(z-Xi)dz+Xi  J  zO(z-X2)dz 

SW  SW 

oo 

+  J  z(XMax  -  2  +  M)0(z  - XMax))dz 

SW 

oo 

-  J  z(XMax-Z  +  ^l)^{z^2)dz 

SW 

oo  oo 

+  0  j  z(l>{z-Xf^g^))dz-a  j  z<}>iz-X2)dz. 

SW  SW 


(5.41) 


94 


Only  two  tenns  need  to  be  reduced.  These  are  the  term  which  is  the  entire 
fourtii  line  and  the  first  term  of  the  sixth  line.  AU  of  the  rest  have  been  reduced 
above. 


SW 


oo 

=  -cr^  J  +  jU  +  cm)i^{u)du 


=  +  iti)y  (S w  -  X^)  -  a^SiSW  -  Xmox)  ■ 


(5.42) 


(T  J  Z^(Z - XMax))dz  =  <7  f  (x^ax  +  M+  <yu)(l>{u)du 

SW 

oo  oo 

=  +  jU )  J  <l>(u)d  M  +  (T^  J  «0(«)dM  (5.43) 

=  <^(XMax  +  M)^(SW-Zai^  )  +  (7^0(SW^^‘^) . 
Adding  all  of  the  reduced  terms  for  this  case  gives 

QpQRBiiSW)  =  -<7^S{^)  -  + )U(T^0(^) + 

+  oxi(xj  +  fi)(x{SW  -Zj)  +  oH^i  +  n)yiSW  -Xi) 

+  Or^dfSW-Zi)  -CT^(Zi  +  nMSW-Xj) 

-  (T^0(SW-Zi)  -  axi^xi  +  pi)a{SW  -z^)  -a^Ziy(SW-Zi) 

+  ozi(z2  +  iJi)oiSW  -Z2)  +  o^ijiSW  -Z2)  (5.44) 

-  ozi (z2  +  n)aiSW  -Z2) -  ©"^(zi  -Z2  -  iJi)y{SW -X2) 

+  cr^5(SW  -  Z2)  -o'^(z2  +  ju)^>(SW-Z2)  -a^(p{SW  -X2) 

+  -O^iXMax  +  iU)r(5W  -  Zn^ax  )  -  W  -  Z^ox  ) 

^^Max  A^)^(SW  —  Z^faj;)  +  O’  0(SW  —  Zai^)  . 


Collecting  like  terms  then  gives  the  final  form  for  BjfSW)  for  the  case  of 
S  W  S  ZAtot  • 
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Bi(SW): 


QpQr 


-a[5(^)  -  -  /i[7(SW)  -  $(5W)j 

+a[5(SW^  xi)  -  ^SW  -  ATj)] 
+/f[y(SW^)  -  0(SW-;ci)] 
+c{5iSW-X2 )  -  0(SW-X2)] 
+/i[7(SW-a:2 )  -  «I>(SW'-a:2 )] 

-ct[«SW  -xu^)-4,(SW-  xjfc)] 

-(j^Ate,  +  /i)[y(SlV  -Xuax)-  flKSiV 


.  (5.45) 


C  Summary 

This  completes  the  derivations  of  the  fonntilas  for  BjCSW).  In  review, 
when  0  <  SW  <  then  BjCSW)  is  given  by  equation  (5.22);  when  Xi  <  SW  <  X2, 
B-iiSW)  is  given  by  equation  (5.30);  when  X2  <  SW  <  Xj^,  Bi(SW)  is  given  by 
equation  (5.39);  and  when  SW  >  Xf^^x'  Bi(SW)  is  given  by  equation  (5.45).  To  get 
B(SW)  we  use  equation  (5.3);  namely, 

BiSW)=Bi{SW) -SWPouriSW) . 

The  equations  for  PQUjiSW)  were  presented  in  Chapter  4. 
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CHAPTER  6  -  MODELS  FOR  ESTIMATING  SAFETY  STOCK 


A.  Introduction 

Hadley  and  Whitin  on  page  165  of  Reference  [3]  state  that  safety  stock  is 
"by  definition,  the  expected  net  inventory  at  the  time  of  arrival  of  a 
procurement."  Their  definition  is,  of  course,  limited  in  applicability  to  only  a 
consumable  item  where  replenishment  of  stock  is  by  procurement.  However,  a 
natural  extension  of  this  definition  to  a  repairable  item  is  possible  when  such  an 
item  is  always  successful  repaired  after  failure.  For  sudi  an  item  there  are  never 
any  procurements  and  safety  stock  can  therefore  be  defined  as  "the  expected  net 
inventory  at  the  time  of  the  arrival  of  a  batch  of  repaired  item." 

Considering  the  pure  prcxnirement  case  first,  the  reorder  point  R  can  be 
defined  as 

R  =  D*PCLT  +  SSp 

where 

D-  expected  quarterly  demand; 

PCLT  =  procurement  lead  time,  quarters; 

SSp  =  procurement  safety  stock; 

and  D  *  PCLT  represents  the  expeded  demand  during  procurement  lead  time. 
Rewriting  this  equation  gives 

SSp  =R-D*  PCLT. 

We  are  concerned  about  representing  SSp  in  terms  of  the  maximum 
inventory  position  for  a  readiness-base  repairable  item  inventory  model.  As  in 
other  chapters  in  this  report  we  will  use  SW  to  represent  the  maximum  inventory 
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position  for  the  readiness-based  sparing  models.  From  Hadley  and  Whitin  [3], 
page  181,  we  know  that 

IPMax-SW  =  R  +  Qp 

for  a  consumable  item  (to  be  called  hereafter  "the  pure  procurement  case"). 
Therefore, 

SSp=SW-D*PCLT-Qp.  (6.1) 

Similarly,  if  we  have  the  pure  repair  case  with  batched  repair  of  quantity 
Qr  then 

SSr=SW-D*  RTAT  -  Qp;  (6.2) 

if  RTAT  is  the  batdi  repair  time. 

If,  on  the  other  hand,  we  have  a  batch  of  Qp  carcasses  being  repaired 
where  each  unit  returns  to  inventory  from  the  repair  depot  as  soon  as  it 
completes  repair  and  there  is  a  time  between  inductions  of  units  into  the  repair 
process,  denoted  by  REP,  we  have  a  special  case  of  Baker's  model  [1]  presented 
in  his  Qiapter  IV  (Non-Instantaneous  Repair  Assessment)  and 

SSp  =  SW  -  D  RTAT  +f  IrEP  -  ^ .  (6.3) 

L  V  2  ;  J  2 

Here  the  average  repair  turnaround  time  is 

RTAT  + 

and  the  average  number  of  units  repaired  over  that  time  is  ^ . 

Equations  (6.1),  (6.2),  and  (6.3)  have  been  confirmed  through  simulation 
runs  made  in  the  process  of  developing  a  generalized  safety  stock  model.  The 
definition  for  the  generalized  safety  stock  was  described  in  Maher  [4]. 
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It  is  "the  expected  net  inventory  at  the  time  of  arrival  of  a  procurement  and/or 
the  arrival  of  a  repaired  carcass  from  the  repair  depot" 

The  question  of  what  form  the  generalized  safety  stock  equation  should 
take  is  not  obvious.  For  a  general  repairable  item  some  repairs  are  successful  and 
some  procurements  are  needed  to  replenishment  the  attritions  due  to 
unsuccessful  repair  or  a  carcass  not  being  returned  for  repair  by  the  customer 
when  he  demands  a  functioning  unit  to  replace  the  failed  one.  Because  both 
procurements  of  new  units  to  replace  old  units  which  could  not  be  repaired  and 
successful  repair  of  other  old  ones  contribute  to  the  ready-for-issue  (RFI) 
inventory,  the  effect  of  the  carcass  return  rate  (CRR)  from  the  customers  and  the 
repair  survival  rate  (RSR)  of  the  repair  process  should  play  an  important  role  in 
die  generation  of  the  general  safety  stock  formula.  A  first  thought  was  that  the 
formula  should  be  a  convex  combination  of  the  pure  procurement  safety  stock 
formula  and  the  pure  repair  safety  stock  formula  where  the  product  CRR  *  RSR 
would  provide  the  weighting  factor;  namely, 

SS  =  (1  -  CRR  *  RSR)SSp  +  CRR  *  RSR  *  SSp ;  (6.4) 

since  0  <  CRR  <  1.0, 0  <  RSR  <  1.0  and  CRR  =  RSR  =  0  for  pure  procurement, 

CRR  =  RSR  =  1.0  for  pure  repair.  Equation  (6.4)  results  in 

SS  =  (1  -  CRR  *  RSR)(SW  -D*PCLT  -  Qp) 

+  CRR  *  RSR(SW  -D*  RTAT  -  Qp) 

=  SW  -  (1-  CRR  *  RSR)D  *  PCLT 

-  CRR  *  RSR  *D  *  RTAT 

-  (1  -  CRR  *  RSR)Qp  -  CRR  *  RSR  *  Qp 

=  SW  -  PPV-lil -CRR  *  RSR)Qp  +  CRR  *  RSR  *Qp]} 

where 

PPV  =  {1  -  CRR*  RSR)D*  PCLT  +  CRR*  RSR  *D*  RTAT .  (6.6) 
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PPV  stands  for  the  Trogram  Problem  Variable/’  the  name  used  by  the  NAVICP 
for  the  expected  demand  during  tiie  aggregate  lead  time  L2,  where 

L2  =  (1  -  CRR  *  RSR)PCLT + CRR  *  RSR  *  RTAT .  (6.7) 

In  the  case  where  we  allow  each  item  entering  repair  to  wait  a  time  REP 
tmtil  its  predecessor  has  passed  through  the  first  workstation,  we  get  the 
modified  form  for  PPV  due  to  Baker  [1]:  namely. 


ZB  =  (1  -  CRR  *  RSR)D*PCLT+CRR  *RSR*D 


RTAT  + 


Qx~'^ 


REP 


(6.8) 


Equation  (6.4)  then  takes  on  the  following  form  based  on  equations  (6.1)  and 
(6.3). 


SS  =  SW-ZB- 


(1-  CRR  *  RSR)Qp  +  CRR  *  RSR 


Qr 

2 


(6.9) 


As  we  will  see  below,  equations  (6.5)  and  (6.9)  are  used  as  the  points  of  departure 
in  the  search  for  other  good  estimates  of  the  safety  stock. 


B.  Simulation  Study 

A  simulation  program  was  written  in  SIMAN IV  that  simvilates  the 
repairable  item  inventory  management  system.  A  description  of  the  SIMAN 
modeling  process  is  provided  by  Reference  [6].  Figure  1  shows  the  details  of  the 
management  process  and  Appendix  A  presents  the  SIMAN  program  used  in  the 
study  and  a  sample  output.  The  simulation  model  assiunes  that  there  is  a  single 
inventory  of  RFI  units,  which  are  ready  for  issue.  Demands  for  RFI  units  arrive 
according  to  a  Poisson  process  and  the  mean  rate  is  D  units  per  quarter.  A 
carcass  accompanies  the  demand  for  an  RFI  unit  with  probability  CRR  and  a 
carcass  is  returned  with  the  demand  with  probability  1.0  -  CRR .  It  is  assumed 
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Figure  1.  Flow  chart  of  the  Repairable  Inventory  Management  Process  [1]. 

that  if  a  carcass  is  not  returned  with  the  demand  then  it  is  damaged  beyond 
repair  and  can  therefore  be  recordable  as  an  attrition.  A  record  of  the  current 
number  of  attritions  is  kept.  When  Qp  attritions  have  been  accumulated  in  this 
record,  a  procurement  action  is  immediately  initiated  for  a  batch  of  Qp  new 
rmits.  These  units  arrive  PCLT  quarters  later.  Returned  carcasses  are 
accumulated  imtil  there  are  Qp  units.  Carcasses  are  then  sent  to  the  depot  in  a 
batch  of  size  Qp.  At  the  depot  the  carcasses  are  queued  and  individually 
evaluated  to  see  if  each  can  be  successfully  repaired.  Each  can  be  successfully 
repaired  witii  probability  RSR.  The  probability  of  failure  is  therefore  1.0-  RSR. 
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Those  tiiat  are  failiires  are  recorded  as  attritions.  When  carcasses  are  inducted  for 
repair,  the  first  one  enters  the  repair  process  immediately.  It  remains  in  repair 
for  RTAT  quarters.  After  it  has  been  in  repair  for  REP  quarters,  the  next  carcass 
is  inspected  and  either  inducted  or  deemed  to  be  an  attrition.  If  it  is  deemed  to 
be  repairable  it  will  be  in  repair  for  RTAT  quarters.  After  it  has  been  in  repair  for 
REP  quarters  the  next  carcass  is  inspected.  As  soon  as  a  carcass  has  been 
repaired,  it  is  returned  as  RFI  and  placed  in  the  inventory.  If  REP  =  0  then  there 
is  no  delay  between  inspections  of  the  individual  carcasses  and  all  repaired 
carcasses  return  to  the  inventory  after  RTAT  quarters. 

The  study  examined  both  the  cases  of  REP  =  0  and  REP  >  0 .  There  were 
70  simulation  runs  with  REP  =  0  and  58  with  REP  >  0 .  Each  run  was  for  1000 
quarters.  Other  parameters  which  were  varied  were  SW,  D,  OO?,  RSR,  PQ.T, 
RTAT,  Qp,  Qp  and  REP  (for  REP  >  0).  In  the  first  part  of  the  study  the  cases  of 
pure  procurement  and  pure  repair  were  simulated.  For  pure  procurement 
equation  (6.1)  was  used  as  the  estimate  in  a  regression  analysis  of  the  simulated 
safety  stock  for  this  case.  As  can  be  seen  in  Figure  2,  the  results  are  excellent. 

Pure  repair  needs  to  be  simulated  under  both  the  conditions  of  REP  =  0 
and  REP  >  0  since  these  situations  do  exist  in  the  Navy.  For  the  case  of  REP  =  0, 
equation  (6.2)  was  used  as  the  estimate  of  the  safety  stock.  For  the  case  of 
REP  >  0,  equation  (6.3)  was  used.  Figures  3  and  4  show  the  respective  results  of 
the  simulations  for  pure  repair. 

The  next  analyses  involved  the  situation  where  REP  =  0  but  where  some 
units  are  attritions  and  others  are  successfully  repaired.  Several  equations  were 
evaluated  for  their  feasibility  as  estimates  of  safety  stock.  Equation  (6.9)  was  used 


102 


Figure  2.  Pure  Procurement. 


as  a  starting  point  for  the  search  for  good  equations  that  would  fit  the  safety 
stock  data  well.  Many  variants  were  tried.  However,  only  those  which  were  as 
good  or  better  than  equation  (6.5)  were  considered  as  potential  alternative 
estimates  of  safety  stock.  The  best  fit  was  obtained  from  equation  (6.10). 

SS  =  SW  -  PPV -  - Q^e-(10-CRJ?*KSR)  _ 

The  other  equations  which  were  tested  and  performed  nearly  as  well 
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were  combinations/ variations  of  equations  (6.5)  and  (6.10).  For  example, 

SS  =  S W  -  PPV  -  Qp  (1.0  -CRR*  RSR)  -  QpRSR ,  (6.11) 

and 

SS  =  SW  -  PPV-Qpe-^^*^^ -QrCRR  *  RSR .  (6.12) 


SS  =  1.0370  +  .98611  *  EQ  (6.2) 
Correlation:  r  =  .99998 


Figure  3.  Pure  Repair  with  REP  =  0. 
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Figure  4.  Ptire  Repair  with  REP  >  0 . 


Figure  5  shows  the  fit  of  equation  (6.5)  to  all  of  the  simulated  safety  stock 
data;  that  is,  all  data  including  pure  procurement  and  pure  repair.  Figttre  6 
shows  the  fit  of  equation  (6.10).  The  regression  equation  for  the  fit  of  equation 
(6.11)  was 

SS  =  -0.6851  +  .98932  *  EQ.(6.11) 
which  has  an  r  =  0.99835. 
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The  regression  equation  for  the  fit  of  equation  (6.12)  was 


SS  =  -1.790  + 1.0082  *EQ.(6.12) 

and  has  r  =  0.99811. 


SS  =  -2.294  +  .99969  *  EQ.  (6.5) 
Correlation:  r  =  .99766 


EQUATION  (6.5) 


Figure  5.  Equation  (6.5)  applied  to  aU  data  for  REP=0. 


06 


SS  =  0.42394  +  .98875  *  EQ.  (6.10) 
Correlation:  r  =  .99957 


EQUATION  (6.10) 

Figure  6.  Equation  (6.10)  applied  to  all  data  for  REP=0. 

Next  we  consider  the  case  where  REP  >  0  and  there  are  both  attritions  and 
repairs.  The  base  case  here  was  given  by  equation  (6.9).  The  regression  results 
are  presented  in  Figure  7.  The  equivalent  of  equation  (6.10)  for  the  case  where 
REP  >0  is: 

SS  =  SW-ZB-  ^  ^3^ 

Figure  8  shows  the  regression  results  when  equation  (6.13)  is  used 
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Figure  7.  Equation  (6.9)  for  all  data  for  REP>0. 


Other  equations,  which  were  examined,  indude 

SS=SW-ZB-Q,(1-  CRR  *RSR)-  ^RSR , 


(6.14) 


and 


SS  =  SW-ZB-Qpe-^^*^^^  -^CRR*RSR.  (6.15) 
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As  can  be  seen,  these  are  similar  to  equations  (6.9)  and  (6.13)  above.  Equation 
(6.14)  regression  results  were 

SS  =  -1.312  +  0.99855  *£Q.(6.14), 

with  r  =  0.99836.  The  regression  results  for  equation  (6.15)  are  shown  below  in 
Figure  9.  This  equation  gives  results  which  are  very  close  to  those  of  equation 
(6.13). 


SS  =  -0.1507  +  .99602  *  EQ.  (6.13) 
Correlation:  r  =  . 99947 


Figure  8.  Equation  (6.13)  for  all  data  for  REP>0. 
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Figure  9.  Equation  (6.15)  for  all  data  for  REP>0. 


Another  equation  provided  an  interesting  different  way  of  making  the 
convex  combination  of  Qp  and  Qp.  It  was 


SS  =  SW-ZB- 


Qp(l  -  CRR  *  RSR)PCLT+^CRR*  RSR{RTAT  +2e^rep) 
(1  -  CRR  *  RSR)PCLT  +  CRR  *  RSR{RTAT  + 


(6.16) 
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This  represents  a  weighting  which  comes  from  the  two  parts  of  the  modified 


form  of  L2  for  tiie  case  of  REP  >  0 . 

L2  =  (1  -  CRR  *  RSR)PCLT + CRR  *  RSRiRTAT + ^^REP) .  (6.17) 

Figure  10  shows  that  the  fit  is  not  as  good  as  the  base  equation  (6.9),  however. 
The  eqxiivalent  form  of  equation  (6.16)  for  the  case  of  REP  =  0  also  gave  a  poorer 
fit  (r  =  0.99728)  than  the  base  equation  (6.5)  for  that  case. 


SS  =  -Z683  +  1.0184  *  EQ.  (6.16) 
Correlation:  r  =  . 99812 


EQUATION  (6.16) 


Figure  10.  Equation  (6.16)  for  all  data  for  REP>0. 


Ill 


C.  Conclusions 


To  place  the  simulation  regression  results  in  the  proper  perspective  it  is 
useful  to  note  that  Professor  Keebom  Kang,  a  simulation  expert  at  the  Naval 
Postgraduate  School,  expressed  amazement  at  how  well  all  of  the  equations 
presented  above  fit  the  simulated  data.  He  noted  that  an  r  >  0.8  is  considered 
very  good  for  such  approximating  equations  in  fitting  simulation  results  for 
complex  processes.  Thus,  the  user  can  decide  on  which  of  these  equations  suits 
his  intuition  and  be  assured  that  he  has  a  good  estimate  for  the  safety  stock  for  a 
given  set  of  parameters. 

It  is  important  to  mention  here  that  the  equations  for  safety  stock 
presented  in  this  diapter  are  applicable  only  for  the  conditions  of  Baker's  model; 
namely,  that  demand  during  the  aggregate  lead  time  is  Poisson  distributed.  In 
the  simulation  model  described  above  the  demand  rate  is  Poisson  and  the 
aggregate  lead  time  is  a  constant  in  PCLT,  RTAT  and  REP.  The  aggregate  lead 
time  L2  may  seem,  however,  to  really  not  be  constant  because  of  the  CRR  and 
RSR  terms  as  equations  (6.7)  and  (6.17)  show.  Forhmately,  Baker  [1]  noted  that 
since  the  individual  parts  are  Poisson  distributed  (i.e.,  demand  during  PCLT  and 
demand  during  RTAT),  the  CRR  and  RSR  terms  are  weights  associated  with  the 
summing  of  two  Poisson  distributed  random  variables  so  that  demand  during 
the  aggregate  lead  time  is  also  Poisson. 

The  conditions  associated  with  the  derivations  in  Chapters  4  and  5  for  the 
assumption  that  the  demand  dxuing  aggregate  lead  time  is  Normally  distributed 
can  be  expected  to  lead  to  different  equations  for  the  safety  stock.  This  will  be 
the  subject  of  a  future  study. 
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CHAPTER  7  -  APPROXIMATE  FORMULAS  FOR  LEAST  COST  Qp  AND  Qp 


A.  Introduction 

This  chapter  presents  the  results  of  a  study  to  determine  approximate 
equations  for  the  values  of  Qp  and  Qp  whidi  minimize  the  average  annual  total 
variable  costs  for  managing  an  inventory  of  an  item.  The  costs  include  the  costs 
to  buy  new  units  (procurement),  the  costs  to  repair  carcasses,  the  costs  to  hold 
items  in  inventory,  and  the  costs  of  backorders.  The  purpose  in  doing  this  was  to 
reduce  the  amoimt  of  computer  time  required  to  determine  the  maximum 
inventory  position  values  which  will  minimize  the  aggregate  mean  supply 
response  time  (MSRT)  for  a  group  of  items  belonging  to  a  weapon  system. 

At  present  the  UICP  model  [5]  uses  the  following  "optimal  values"  for  Qp 
and  Qp. 


where, 

D  =  Quarterly  Demand  Rate; 

A  =  Procurement  Contract  Costs; 

A2  =Repair  Contract  Costs; 

C  =Unit  Purchase  Cost; 

C2  =  Unit  Repair  Cost; 

I = Holding  Cost  Rate  =  0.21; 

G  =  CRR*RSR*D; 

CRR  =  Probability  of  a  Carcass  Being  Returned  by  a  Customer; 
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RSR  =  Probability  of  a  Carcass  being  Repaired. 

These  equations  were  derived  from  conjectured  UICP  average  annual  total 
variable  costs  equations,  one  for  managing  newly  procured  units  and  one  for 
managing  repairing  of  carcasses  and  the  inventory  of  successful  repaired  units 
The  new  repairable  model  has  a  cost  equation  which  combines  the  two 
parts  conjectured  by  the  UICP  modelers  since  the  repaired  units  and  newly 
procured  units  are  kept  in  the  same  inventory  as  equally  valid  ready-for-issue 
(RFI)  units.  This  equation  can  be  written  as 


JVC  =  +  +  ic^eOH  +  mSW) ,  {72) 

s2p  Qr 

where  EOHrepresents  the  expected  on-hand  inventory.  Since  the  expected  net 
inventory,  E(NI),  is  the  difference  between  the  expected  on-hand  and  the 
expected  number  of  backorders,  B(SW),  we  can  write  EOH  -  E{NI)  +  B{SW), 
which,  upon  substitution  of  the  E(NI)  formula  of  Baker  [1],  gives 


EOH  =  SW-ZB  +  e{SW)  . 


Here, 


SW  =  Maximum  Inventory  Position, 

ZB  =  Baker’s  Extension  for  PPV  (given  by  equation  (6.8)). 

In  equation  (7.2)  the  average  annual  number  of  repair  contracts  is  given  as 


4  *  CRR *  D 
Or 


.  The  number  of  accumulated  carcasses 


per  year  is  the  numerator 


since  CRR  represents  the  probability  of  a  carcass  being  returned  per  replacement 
unit  demanded.  The  denominator  is  Qj^  since  that  is  the  number  of  carcasses 
inducted  as  a  batch  (Qj^  does  not  represent  the  number  of  carcasses  repaired). 


114 


The  expected  aimual  backorder  cost  term  in  equation  (7.2)  is  given  by 
AB(SW).  Formulas  for  computing  B{SW)  were  derived  in  Qiapters  3  and  5  of 
this  report.  The  shortage  cost  per  year  per  unit  is  represented  by  A.  In  the  UICP 
model  the  shortage  cost  is  per  requisition  and  there  is  also  an  essentiality  term 
added  into  the  backorder  costs.  Thus,  if  such  a  form  is  of  interest  then  we  can 

write  the  equivalence  as  A  =  ^  '%1CP  where  E  r^resents  the  niimber  for 


essentiality  (currently  0.5  for  Mechanicsbuig)  and  F  represents  the  average 
requisition  size.  In  the  new  model  F  =  1.0  is  assumed  (a  review  of  repairable 
item  data  shows  this  to  be  approximately  true). 

Finally,  ttie  holding  cost  rate  in  equation  (7.2)  includes  the  average  unit 
cost,  called  C3.  which  is  a  convex  combination  of  the  xmit  purchase  cost  and  the 
unit  repair  cost. 


C3  = 


C2. 


(7.3) 


C.  Computer  Runs 

A  computer  program  was  written  in  FORTRAN  77  to  determine  tiie 
optimal  values  for  Qp  and  Qp  which  will  minimize  equation  (7.2)  for  a  wide 
range  of  parameter  values.  These  included  D,  CRR,  RSR,  PCLT,  RTAT,  REP, 
C,  C2,  A,  A2,  and  A. .  The  program  was  designed  to  run  for  a  fixed  set  of  these 
parameters  over  a  range  of  SW  values  from  0  to  150  or  200  in  increments  of  10. 
Sixty-five  such  sets  were  obtained. 

The  approach  taken  in  the  program  is  to  incrementally  increase  Qp 


starting  at  1  and  then  incrementally  increase  Qp  to  find  the  optimal  Qp,  Qp,  for 
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each  Qj^  value.  As  long  as  TVCiQp,Qj^)  is  decreasing  the  value  of  Qp  is 
increased.  This  search  process  continues  \mtil  the  next  value  of  Qp  shows  an 

increase  in  TVCiQp,Qp ).  Optimal  Qp  is  then  the  previous  value  of  Qp . 

The  result  of  these  computer  runs  were  tiien  recorded  in  a 
STATTSTTCA  4.1  [8]  file.  Files  were  developed  for  a  range  of  values  of  the 
product  CRR  *  RSR  and  each  of  the  other  parameters. 

Only  one  parameter  was  allowed  to  vary  at  a  time.  Plots  of  the  Qp  and 
Qp  values  as  a  function  of  a  parameter  were  examined  to  see  what  function  of 
that  parameter  would  provide  a  good  regression  fit.  These  were  often  quite 
complex  nonlinear  functions.  When  the  functions  had  all  been  determined,  they 
were  assumed  to  be  combinable  into  a  sum  of  separate  functions  ( cross  products 
and  other  combinations  were  examined  prior  to  deriding  on  this  approach) 
which  was  then  applied  to  the  complete  set  of  data  to  see  how  good  the 
aggregate  fit  was.  The  results  were  surprising  good. 

In  the  process  of  conducting  the  analyses  it  was  discovered  that  when 
SW —ZB  <  0.0  that  the  optimal  values  of  Qp  and  Qp  remained  unchanged  as  SW 
was  increased.  Figure  11  shows  ttie  t3q)ical  behavior  of  optimal  Qp  and  Qp  as  a 
function  of  SW -ZB. 

Mathematical  analyses  of  TVC  to  see  what  formulas  could  be  derived  for 
unit  changes  in  Qp  and  Qp  were  conducted.  The  first  analysis  was  to  find  the 
relation  for  determining  optimal  Qp .  Using  the  technique  of  finite  differences,  we 
seek  the  largest  value  of  Qp  for  which 

ATVCiQp)=  TVCiQp)-TVC{Qp  - 1) <  0. 
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Figure  11.  Optimal  Qp  and  as  a  function  of  SW  -  ZB. 
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4A(D-G)  4AyCRR*D 

ATVQQp )  = 


+ia 


SW  -  D(1  -CRR  *  RSR)PCLT 
-D*CRR*RSR 

Qp  +Qr 


RTAT  +  ^^REP 

2  j 


+  msw,Qp)- 


■  +  B(SW,Qp) 
4A(D-G)  4A2CRR*D 


Qp-1 


Qr 


-IC. 


S  W  -  D(1  -CRR*  RSR  )PCLT 
-D*CRR*  RSr(  RTAT  +  ^^REP 

I  2 


(Qp  ~1}  +  Qp  ~2 

•  -  \  —  +  B(SW,Qp  -1) 

-XB[SW,Qp  -1). 


(7.4) 

After  collecting  tenns  and  canceling  common  terms  in  equation  (7.4)  we  have 
ArVC(Qp )  =  -  -  +  (IC3  +  A)AB(S  W,  Qp  )  S  0,  (7.5) 

where 

AB{SW,Qp)  =  BiSW,Qp)-B{SW,Qp  -1). 

From  inequality  (7.5)  we  can  write  the  optimality  condition  as 

®'’(®P‘®^2(JC3  +  A)AB(SW,Qp)-7C3  ' 

No  attempt  was  made  to  use  the  formulas  for  B{SW)  from  Qiapter  3 
because  of  their  complexity.  Instead,  a  detailed  investigation  of  the  computer 
results  was  conducted.  It  showed  that  for  SW  <  ZB  that  AB{SW,Qp)  =  0.5  in  the 
neighborhood  of  optimal  Qp.  Therefore,  inequality  (7.6)  reduces  to 
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(7.7) 


Qp(Qp  -1)^ 


8A(D-G) 

X 


Thus,  we  want  the  largest  value  of  Qp  for  which  inequality  (7.7)  is  satisfied.  A 
check  of  the  computer  results  also  showed  that  the  optimal  Qp  from  the 
computer  search  procedure  was  that  value  which  satisfied  inequality  (7.7)  when 
SW  <  ZB. 


To  find  an  inequality  for  determining  optimal  Qpthe  approach  used  was 
similar  to  that  for  Qp.  Now  we  want  the  largest  value  of  Qp  for  which 

atvc{Qr) = ryc(Qp)  -TvciQ^  - 1)  <  o . 

In  this  case 


A  A  (^T?  jf?  ^  TC^ 

AiraQs)  =  ^  —f{l+G»REP}+  {ICi  +  ^)AB(SW,Qs) .  (7.8) 

The  computer  results  in  the  neighborhood  of  optimal  Qp  were  not  a  constant  like 
was  the  case  for  Qp.  Values  were  recorded  in  a  STATISTICA.  4.1  file  [8]  for  a 
variety  of  parameters  and  their  values.  Regression  analyses  discovered  that 

AB(SW,Qp)  =  |(1 +G  »  REP).  (7.9) 


Figure  12  shows  the  impressive  fit  of  the  data  to  equation  (7.9).  As  a 
consequence,  the  optimality  condition  for  Qp  when  SW  <  ZB  is  that  it  is  the 
largest  value  of  Qp  for  which 


Qr(Qr 


8A2CRR  *  D 
X(1+G*REP)  ' 


(7.10) 


Inequality  (7.10)  was  confirmed  by  the  computer  program  results. 


119 


When  SW  >  ZB  it  is  easy  to  see  that  optimal  Qp  and  Qp  are 
approximately  linear  in  the  difference  SW  -  ZB.  However,  an  attempt  to 
develop  some  relationship  between  AB(SW,Qp)  or  AB{SW,Qp)  and  optimal  Qp 
or  Qp  was  unsuccessful.  The  observed  values  of  AB{SW,Qp)  and  AB{SW,Qp) 
become  erratic  and  quite  small  in  value  as  SW  increases.  Ottier  terms  of  the 


AB(SW,  QR)  =  0.00092  +  .99983  *  EQ.  (7.9) 
Correlation:  r  =  . 99999 


AB(SW,QR) 


EQUATION  (7.9) 


Figtu-e  12.  The  Regression  Fit  of  Equation  (7.9)  to  AB(SW,  Qp)  Data. 
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TVC  equation  also  dominate  the  detennination  of  optimal  Qp  and  Qp  as  SW 
increases.  As  a  consequence,  an  investigation  of  the  relationships  between 
optimal  Qp  and  Qp  and  each  of  the  parameters  (except  I)  of  TVC  was  begun. 

As  mentioned  earlier,  the  process  consisted  of  varying  the  value  of  a 
parameter  while  all  of  the  others  remained  fixed.  From  plots  of  the  values  of 
optimal  Qp  or  Qp  for  a  range  of  a  parameter's  value  functions  of  the  parameters 
were  suggested  and  fits  attempted.  Then  improvements  were  made  and  tested. 
This  iterative  process  eventually  closed  in  on  the  best  fitting  function. 

Figure  11  shows  optimal  Qp  and  Qp  to  be  a  linear  function  of  SW  -  ZB . 
We  know  from  the  equation  for  ZB  given  by  equation  (6.8)  that  it  contains  D, 
CRR,  RSR,  PCLT,  RTAT,  REP,  and  Qp.  However,  there  are  more  fuiictions  of 
these  and  other  parameters  which  affect  the  slopes  and  intercepts  of  those 
straight  lines  when  SW  >  ZB. 

To  illustrate  the  process,  we  will  use  the  analysis  of  Cs's  effect  on  optimal 
Qp.  We  first  define 

DELQP^Q;-Q;iSW  =  0), 

where  Qp(SW  =  0)  satisfies  equation  (7.7)  when  SW  =  0. 

Figure  13  shows  curves  of  DELQP  as  a  function  of  SW  -  ZBfor  the  range 
of  C3  values  from  $720  (the  lower  curve)  to  $10,800  (the  upper  curve).  Since  we 

already  know  how  to  compute  Qp(SW  =  0)  we  focus  on  SW  >  ZB. 

Figtire  13  shows  the  plots  of  DELQP  are  fairly  linear  in  S  W  -  ZB  but  the 
slopes  and  intercepts  change  as  C3  increases.  Figtjres  14  and  15  show  the  slopes 
and  intercepts  for  the  extremes  of  the  C3  values  resulting  from  regression  fits. 
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For  each  of  tiie  curves  of  Figure  13  the  slopes  and  intercepts  were  determined  as 
was  done  in  Figures  14  and  15.  These  slope  and  intercept  values  were  then 


Figure  13.  Optimal  Values  of  DELQP  for  a  Range  of  C3  Values. 


plotted  as  a  function  of  C3  as  shown  in  Figures  16  and  17.  Using  a  lot  of  guess 
work,  various  functions  were  conjectured  and  tried  until  something  reasonably 
close  to  eadi  of  the  curves  in  Figures  16  and  17  was  foxmd.  In  Figure  16  QPSLP 
stands  for  the  slope  of  the  Qp  curve  for  a  certain  C3  value.  In  Figure  17  QPESJT 
stands  for  the  intercept  of  the  Qp  curve  at  SW  =  ZB . 
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For  QPSLP.  the  function 


f: 


\022 


Viooo; 


was  fotmd  to  be  dose.  A  regression  fit 


with  QPSLP  suggested  that  a  coeffident  of  0.554  and  an  intercept  constant  of 
-0.031  was  needed.  Thus,  the  equation  which  will  be  used  as  an  estimate  for 
QPSLP  is 


QPSLPiC^)^  -0.031+  0.554|^' 


r  \022 

• 

1000 J 


(7.11) 


DELQP  =  -4.054  +  .90894  *  (SW  -  ZB) 
Correlation:  r  =  . 99919 


DELQP 


Figure  15.  Regression  Fit  of  DELQP  with  SW  -ZB  for  C3  =  $10, 800  when 

SW>ZB, 

Figure  18  shows  the  regression  fit  between  QPSLP  and  equation  (7.11).  The 
value  of  the  correlation  coefficient  r  is  0.98092  shows  tiie  fit  is  good.  Figure  19 
shows  how  the  two  of  them  compare  when  they  are  plotted  as  functions  of  C3 . 

The  intercept  part,  QPINT,  was  more  difficult  to  develop  a  function  for. 
The  equation  which  was  finally  chosen  is 
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QPINTiC^)  =  -5.86  +  0.18 


2 


(7.12) 


The  regression  fit  gave  r  =  0.94565  which  is  not  as  good  as  equation  (7.11)  for 
QPSLP.  The  comparative  plots  of  equation  (7.12)  and  QPINT  as  a  function  of  C3 
are  shown  in  Figure  20.  Qearly,  equation  (7.12)  provides  only  approximate 
results. 


Figure  16.  Sopes  of  C3  Curves  for  Optimal  Qp  when  SW  >  ZB. 
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Figtire  17.  Intercepts  of  C3  Cvtrves  for  Optimal  Qp  when  SW  >  ZB. 

The  other  approximate  formulas  for  the  slopes  and  intercepts  for  Qp  for 
the  other  parameters  are  given  in  Table  1.  Table  2  provides  comparable 
information  for  Qp. 
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Figiire  18.  Regression  Results  from  Fitting  Equation  (7.11)  to  QPSLP  Data. 
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Figure  19.  Comparison  of  Equation  (7.11)  to  QPSLP  Data  for  all  C3  Values. 


For  both  Qp  and  Qp  the  effects  of  PCLT  and  RTAT  were  completely 
described  by  ZB  so  fliere  are  no  formulas  in  the  tables  for  these  parameters.  It  is 
important  to  note  that  the  formulas  involving  REP  did  include  REP  =  0,  which 
corresponds  to  batch  repair,  in  their  development. 

In  the  combining  of  the  various  effects  the  summing  all  of  the  constant 
terms  gave  values  which  was  far  from  what  should  be  the  values  to  give  a 
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reasonable  slope  and  intercept  for  Qp  and  Qp  for  any  set  of  values  of  parameters 
and  a  range  of  SW  values  from  ZB  to  150.  Therefore,  those  constants  were 
ignored  and  constants  were  determined  by  selecting  several  sets  of  data  and 
comparing  their  known  slopes  and  intercepts  with  those  being  provided  by  the 
sum  of  effects  without  constants  included.  An  average  value  was  then  selected 
for  the  summed  effects’  constants.  For  Qp  the  combined  effects’  slope  constant 
was  -2.392  and  the  intercept  constant  was  -7.1. 
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Figure  20.  Comparison  of  Equation  (7.12)  to  QPINT  Data  for  all  C3  Values. 
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Table  1.  Slope  and  Intercept  Equations  for  Estimating  Optimal  Qp. 


Qp  SLOPE  &  INTERCEPT 

Correlation 

QPSLPiD,  CRR  *RSR)  =  -1.04  +  1.67D°°^ 

+ 1^- 1 .602 + ^  jcRR  *  RSR 

r=  0.99268 

QPINTiD,  CRR*  RSR)  =  -0.297  +  D(0.19  -  0.008D) 

- D(1.76  -0.08D)iCRR*  RSR  -  0.060)^ 

r=  0.78815 

QPSLPiX)  =  0.01157 + 

r=  0.98933 

QPINTW  =  0.661+5.208»lnfY^l 

r=  0.95560 

QPSLPiREP)  =  0.6567  -  0.2927REP®^^ 

r=  0.99621 

QPINTiREP)  =-3.508  +7288REP^-^ 

r=  0.99020 

QPSLPiA)  =  0.469  +  0.00003 A 

r=  0.98735 

QPimiA)  =0.622  -  0.00026A 

r= -0.5735 

QPSLPiA2)  =  0.5257  +  O.OOOOIA2 

r=  0.81992 

QPINTiA2)  =  3.2027  -  0.00356A2 

r= -0.9805 

130 


Table  2.  Slope  and  Intercept  Equations  for  Estimating  Optimal  . 


Qr  slope  k  intercept 

Correlation 

QRSLPiC^)  =  0.31+0.45 

aoooj 

r=  0.99152 

QRINTiCs)  =-0.0656- 

r=  0.96689 

^55.1  +  0.0195C3j 

QRSLPiD,CRR  *  RSR)  =  0.4  +  0.36(CRR  *  RSRf^^  +  0.045Vd 

r=  0.95497 

QRINTiD.CRR  *  RSR)  =-4.52  +  d[-0.875 + 1.8(CRR  *  RSR  -  0.35)^] 

r=  0.99776 

QRSLP{X)=  0.756+^^ 

A 

r=  0.99912 

z'  i  _  7nn 

QRINTiX)  =  -12.22  -  0.066|^ J 

r=  0.89309 

QRSLPiREP)  =  0.6337  +  0.3409REP®-^ 

r= -0.9933 

QRINriREP)=  -3.508+  7.288PEP^^ 

r=  0.99020 

A 

QPSLP(A)  =  0.746  +0.213e  looo 

r=  0.98816 

QRim{A)=  -10.05-  0.00196 A 

r= -0.9794 

QRSLP{A2  )  =  0.80614  -  0.00006 A2 

r= -0.9898 

QPINT(A2  )  =-15.45  +  O.OO253A2 

r=  0.99356 
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The  formula  for  estimating  optimal  Qp  when  SW  >  ZB  is 
QPEST  =  Qp  {SW  =  0) + XQPJNT  +  [XQPSLP](SW  -ZB),  (7.13) 

where  the  Qp  value  used  in  ZB  when  REP  >  0  is  the  value  of  QREST  obtained 
from  equation  (7.15).  The  derivation  of  equation  (7.15)  begins  witii 


QREST  =  Qp  (SW  =  0)  +  QRINT  +[X  QRSLP]{SW  -  ZB(QREST)} 

=  Q*p(SW  =  0)+'^QRINT+[£QRSLP^W 

-  [5;  QRSLP]D(1  -CRR  *  RSR)PCLT 

-  [  J^QBSLPjD  *  CRR  *RSR  *RTAT  (7.14) 

+  [£QRSLPY:RR  *  RSR^ 

-  E  Q^SLPJd  *  CRR  *  RSR  ^  QREST . 


Rewriting  equation  (7.14)  by  moving  QREST  to  the  left-hand  side  of  the 
equation. 


QR£Sr  +  [5;QRSLP]D  *CRR  *RSR-^QREST 

=  (Tr{sw  =  o)  +  J^qrint  +  [J^qrslp}w 

-  [XQ^SLP]D(1-CRR*  RSR)PCLr 

-  [£QRSLP]D  *  CRR  *  RSR  *  RTAT 

+  [XQ^SLPJcRR  *  RSR^ 

=  Qr  (SW  =  0)  +  J^QRim  +  QRSLPjjsW  -  PPV  +  CRR  *  RSR —1 . 


Finally,  the  formula  for  QREST  is 
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QREST= 


Ql (SW  =  0)  +  £  QRINT  +[£QRSLP]|sW  -  PPV  +CRR  *  RSR^j 

^  r^r  T\1  -f-x.  1X/X-IX  K,1£P 


l+[2^QRSIi>]D*CRR*  RSR 


(7.15) 

The  combined  effects’  constants  for  ttie  summed  QRSLP  and  QRINT  effects 
were  0.142  and  6.23,  respectively. 

Figures  21  and  22  show  the  regression  results  when  the  computed  and 
estimated  values  of  optimal  Qp  and  Qp  for  the  65  sets  of  data  are  compared. 


C.  Conclusions 

The  formulas  developed  in  this  chapter  are  much  more  complex  that  those 
currently  by  the  UICP  model.  However,  they  do  give  the  minimum  average 
annual  total  variable  costs  to  manage  an  inventory  of  a  repairable  item  over  a 
broad  range  of  maximum  inventory  position  values.  They  have  a  constant  value 
when  SW  ^  ZB  and  increase  approximately  linearly  with  SW  when  SW  >  ZB. 
The  UICP  formulas  (equation(7.1))  do  not  change  with  SW.  Indeed,  they  are 
only  approximately  optimal  for  one  value  of  the  maximum  inventory  for  any  set 
of  parameters.  Figure  23  illustrates  this  fact  with  the  curves  of  comparative  TVC 
values  for  a  set  of  parameter  data. 

These  formulas  were  developed  to  elimmate  the  requirement  to  conduct  a 
time-consuming  search  for  optimal  Qp  and  for  each  value  of  SW  when 
conducting  a  marginal  analysis  while  trying  to  find  a  set  of  SVV^,  z  =  1,  n  which 
will  minimize  the  aggregate  MSRT  of  a  set  of  n  repairable  items. 

Unfortunately,  cost  minimizing  and  MSRT  minimizing  are  two  different 
objectives.  As  one  expects,  increasing  SW  for  a  fixed  set  of  Qp  and  Qp  values 
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will  reduce  the  MSRT  for  an  item.  However,  changing  Qp  and  Qp  so  that  they 
minimize  TVC  as  SW  changes  causes  the  MSRT  to  increase  as  SW  increases 
beyond  SW  =  ZB.  Figures  24  and  25  show  how TVC-optimal  Qp  and  Qp  change 
with  SW  and  how  MSRT  changes  with  SW  when  these  TVC-optimal  Qp  and  Qp 
are  used. 


Figure  21.  Regression  Results  from  the  Comparison  of  the  Approximate  Estimate 
of  Optimal  Qp  and  the  Computed  Value  of  Optimal  Qp. 


A  companion  report  will  show  that  attempting  to  use  these  TVC-optimal 
Qp  and  Qp  values  in  the  marginal  analysis  will  result  in  a  'less  than  optimal" 
aggregate  MSRT  because  of  the  behavior  shown  in  Figure  25.  If,  instead,  Qp  and 
Qp  are  fixed  prior  to  tihie  marginal  analysis,  then  NCRT  continues  to  decrease  as 
SW  increases  and  an  optimal  MSRT  value  can  be  obtained. 


Figure  22.  Regression  Results  from  the  Comparison  of  the  Approximate  Estimate 
of  Optimal and  the  Computed  Value  of  Optimal  Qp. 
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Figure  23.  A  Comparison  of  Optimal  TVC  and  TVC  having  the  UICP  Values  of 

Qp  and  Qp. 
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Figure  24.  Min  TVC  -Optimal  Qp  and  Qp  for  the  Same  Data  as  Figure  23. 
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Figure  25.  Mmitnum  MSRT  for  Optimal  Qp  and  Qp  of  Figure  23. 
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CHAPTER  8  -  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


A.  Summary 

Chapter  1  reviewed  the  history  of  the  development  of  a  new  repairable 
inventory  management  model  for  the  replenishment  of  stodc,  including  the 
procurement  of  new  units  and  repair  of  damaged  carcasses.  The  goal  of  this 
effort  is  an  analytical  model  of  the  Navy’s  repairable  item  inventory  management 
process  which  will  allow  the  determination  of  stodcage  depths  which  will  meet 
readiness  goals.  This  report  presents  the  next  logical  steps  in  that  development 
process. 

Chapters  2  and  4  presented  the  derivations  of  formulas  for  the  probability 
of  being  out  of  stock  at  any  instant  of  time.  Pout  (SW)/  for  the  cases  of  demand 
during  aggregate  lead  time  being  Poisson  and  Normally  distributed.  These 
derivations  provided  formulas  for  all  possible  values  of  the  maximum  inventory 
position,  SW .  Th^  assumed  that  values  for  the  order  and  repair  quantities,  Qp 
and  Qp,  are  given. 

Chapters  3  and  5  presented  the  derivations  of  formulas  for  the  expected 
number  of  backorders  at  any  instant  of  time,  B(SW),  for  the  cases  of  demand 
during  aggregate  lead  time  being  Poisson  and  Normally  distributed.  As  with  the 
PouriSW)  formulas,  the  derivations  provided  formxalas  for  all  possible  values  of 
S  W  and  assumed  that  the  values  of  Qp  and  Qp  are  given. 

Chapter  6  presented  the  results  of  a  simulation  study  to  determine  an 
approximate  formula  for  the  safety  stock  for  repairable  items  for  the  case  where 
demand  during  the  aggregate  lead  time  is  Poisson  distributed.  Safety  stock  is 
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defined  as  the  expected  net  inventory  at  the  time  of  the  arrival  of  a  procurement 
order  and/ or  the  return  of  a  repaired  carcass  from  a  depot.  Four  formulas  were 
developed;  two  for  the  case  of  batched  repair  (no  delays  between  inductions  of 
units  in  the  batch),  and  two  for  the  case  of  a  delay  between  units  being  inducted. 
These  formulas  are  also  dependent  on  knowing  the  values  of  SW,  Qp  and  Qp. 

Chapter  7  presents  the  approximate  least  cost  formulas  for  Qp  and  Qp  for 
the  case  where  demand  during  the  aggregate  lead  time  is  Poisson  distributed. 
The  formulas  incorporate  the  effects  of  all  model  parameters.  They  are  much 
more  complex  that  tire  UICP  formulas  [5]  and  do  increase  linearly  with  SW 
when  SW  is  greater  than  the  expected  demand  during  the  aggregate  lead  time. 
The  purpose  of  developing  these  formulas  was  to  eliminate  the  time  which 
would  be  required  in  determining  their  values  dining  each  step  of  a  marginal 
analysis  to  find  a  set  of  optimal  SW 's  for  a  group  of  items  which  minimize  the 
aggregate  MSRT  value. 

B.  Conclusions 

Formulas  for  Pqijj{SW)  and  B{SW)  based  on  Baker’s  model  [1]  and  an 
extension  of  it  to  the  case  of  Normally  distributed  demand  during  aggregate  lead 
time  have  been  developed  which  should  be  easy  to  use  in  the  Navy’s  Uniform 
Inventory  Control  Program  (UICP)  [5].  They  are  essential  to  the  implementation 
of  a  readiness-based  replenishment  model.  However,  the  formulas  do  depend 
on  knowing  the  maximum  inventory  position  SW,  the  order  quantity  Qp,  and 
the  repair  quantity  Qp . 

Four  approximate  formulas  for  safety  stock  have  been  developed  which 
are  easy  to  use  and  intuitively  appealing.  However,  they  also  require  knowing 
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SW,  Qp,  and  Qp,  and  d\e  foraiulas  of  Chapter  6  are  only  valid  for  the  Baker 
model. 

The  least  cost  formulas  for  the  order  and  repair  quantities  are  complex  but 
do  eliminate  the  need  for  an  elaborate  search  routine  to  determine  their  values, 
for  a  given  set  of  repairable  item  parameters.  The  major  drawback  to  their  use  in 
a  readiness-based  model  is  that,  because  they  increase  with  SW,  they  tend  to 
cancel  out  the  benefit  (in  the  sense  of  reducing  mean  supply  response  time 
(MSRT))  which  can  be  realized  from  using  a  fixed  pair  of  Qp  and  Qp  values. 

C.  Recommendations 

The  formulas  for  Pqup(SW)  and  B(SW)  are  ready  to  be  used  by  the  UICP 
[5].  The  problem  which  remains  to  be  solved  is  how  to  determine  the  actual 
maximum  inventory  position  (IP)  for  repairable  item.  A  consumable  item’s 
maximum  IP  value  can  be  determined  for  the  sum  of  its  reorder  point  and  order 
quantity.  This  approach  cannot  be  used  for  a  repairable  because  it  is  impossible 
to  determine  a  reorder  point  and  it  is  not  dear  how  the  effects  of  Qp  and  Qp 
should  be  handled.  Fortunately,  an  approximate  value  for  SW  can  be  found  by 
using  one  of  the  formulas  for  safety  stock. 

Unfortunately,  approximate  formulas  for  safety  stock  for  the  case  where 
the  demand  during  the  aggregate  lead  time  is  Normally  distributed  have  not 
been  derived  yet.  This  case  creates  a  problem  for  simulation  modeling  because 
the  probability  distribution  for  the  time  between  demands  is  not  known.  One 
possible  approach  is  to  generate  the  quantity  of  demands  for  a  quarter  from  the 
Normal  distribution  and  then  assume  the  arrival  times  are  equally  spaced  over 
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the  quarter.  Hopefully,  if  enough  quarters  are  simulated  the  result  will  be  a 
good  approximation  to  the  assumed  demand  distribution. 

The  least  cost  formulas  for  Qp  and  Qp  are  not  useful  for  a  readiness-based 
replenishment  model.  As  Figure  25  at  the  end  of  Chapter  7  shows,  MSRT 
increases  as  the  value  of  SW  increases  once  SW  exceeds  the  expected  demand 
during  the  aggregate  lead  time,  ZB.  Fixing  the  Qp  and  Qp  values  at  some  value 
allows  the  value  of  MSRT  to  continue  to  decrease  as  SW  increases.  A  study  of 
the  effects  of  several  different  fixed  values  of  Qp  and  Qp  is  therefore  needed. 

The  results  of  a  preliminary  study  of  candidate  Qp  and  Qp  values  are  presented 
in  a  companion  report. 
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APPENDIX  A  -  SIMAN IV  [6]  PROGRAM  FOR  SAFETY  STOCK 


A.  MODEL  FRAME 

BEGIN,  Repairable  Model  No.4; 

/ 

next  CREATE:  EX(1); 

COUNT:  Demanded; 

ASSIGN:  Inventoiy_position=Inventory_position-l ; 
ASSIGN:  Net_inventory=Net_inventory-l; 

BRANCH,  1:  WITH,  CRR,  carcass: 

ELSE,  no  carcass; 

/ 

no  carcass  ASSIGN:  OrderQ=OrdeiQ  + 1; 

BRANCH,  1:IF,  OrderQ.GE.Q,  order: 

ELSE,  continue; 

/ 

order  ASSIGN:OrderQ=0; 

TALLY:IP  at  order,Inventoiy_position; 
ASSIGN:Inventory_position=Inventoiy_position+Q; 
DELAY:  PCLT; 

TALLY:  Net  inventory  at  receipt,  Net_inventory; 
COUNT:  Ordered,Q; 

ASSIGN:  Net_inventoiy=Net_inventory+Q:DISPOSE; 
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/ 


carcass  QUEUE,  Carcass; 

GROUP:  QR; 

ASSIGN:  lnduct=0; 

ASSIGN:  Repaired=0; 

TALLY:  IP  at  order,Inventory_position; 

ASSIGN:  Inventoiy_position=Inventory_position+QR; 
SPLIT; 

ASSIGN:  Induct=Induct+l; 

BRANCH,  1:  WITH,  RSR,  success: 

ELSE,  bad  carcass; 

/ 

success  ASSIGN:  Repaired=Repaired+l; 

DELAY:  (Repaired-1)*DELTA; 

BRANCH,  1:IF,  DELTA.LE.0.0,  repair2: 

ELSE,  go  on; 
goon  DELAY: RTAT; 

TALLY:  Net  inventory  at  receipt,  Net_inventory; 
COUNT:  Repair; 

ASSIGN:  Net_inventory=Net_mventory+l:DISPOSE; 

/ 

repair2  BRANCH,  1:IF,  Induct.EQ.QR,  repairS: 
ELSE,continue; 
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repairs  ASSIGN:  TempNR=Repaired; 

DELAY:  RTAT; 

TALLY:  Net  inventory  at  receipt,  Net_inventory; 

COUNT:  Repair,  TempNR; 

TALLY:  Repair  batch  size,  TempNR; 

ASSIGN:  Net_inventoiy=Net_inventory+TempNR:DISPOSE; 

bad  carcass  DELAY:  Repaired*DELTA; 

ASSIGN:  Inventory_position=Inventoiy_position-l; 
BRANCH,  1:  IF,DELTA.EQ.O.O.AND.Induct.EQ.QR,  check: 
ELSE,  no  carcass; 

/ 

check  BRANCH,  1:IF,  Repaired.GT.O,clone: 

ELSE,  no  carcassl; 

no  carcassl  ASSIGN:  TempNR=0; 

TALLY:  Repair  batch  size,TempNR:NEXT(no  carcass); 

/ 

done  DUPLICATE:  l,repair3:NEXT(no  carcass); 

/ 

continue  DELAY:0:  DISPOSE; 

/ 

END; 
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B.  EXPERIMENTAL  FRAME 


BEGIN; 

PROJECT,  Repairable  Model  No.4; 

ATTRIBUTES:  TempNR; 

QUEUES:  Carcass:  Nocarcass; 

PARAMETERS:  1,.05967; 

VARIABLES:  Demand,  16.76: 

PCLT,  6.07: 

RTAT,  1.28: 

DELTA,  0.: 

CRR,  .9764: 

RSR,  .85: 

0,6: 

QR,  16: 

OrderQ,  0: 

RepairQ,  0: 

Induct,  0: 

Repaired,  0: 

Num,  0: 

Inventory_position,  59: 

Net_inventory,  59; 

TALLIES:  Net  inventory  at  receipt:  Repair  batch  size:  IP  at  order; 

DSTATS:  Net_inventory,  Expected  Net  Inventory:Inventory_position,  Exp  IP; 
COUNTERS:  Demanded„No:Ordered„No:Repair„No; 
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REPLICATE„4050„,50; 

END; 


C  EXAMPLE  OUTPUT  LISTING 


SIMAN IV  -  License  #9050352 
Naval  Postgraduate  School 

Summary  for  Replication  1  of  1 

Project:  Repairable  Model  No.4  Run  execution  date :  8/9/ 1999 

Replication  ended  at  time  :  1050.0 
Statistics  were  cleared  at  time:  50.0 
Statistics  accumulated  for  time:  1000.0 

TALLY  VARIABLES 


Identifier 

Aver.  Var. 

Min. 

Max.  No.  of  Obs. 

Net  inventory  at  recei  8.7921  .76333 

-22.00 

33.000 

1462 

Repair  batch 

size  13.670  .10017 

8.0000 

16.000 

1007 

IP  at  order 

44.120  .13091 

38.000 

53.000 

1461 

DISCRETE-CHANGE  VARIABLES 

Identifier 

Aver.  Var. 

Min. 

Max.  Final  Value 

Expected  Net  Inventory  14.876  .48251 

-22.00 

39.000 

28.000 

Exp  IP 

49.034  .10005 

39.000 

59.000 

55.000 
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COUNTERS 


Identifier 


Demanded 

Ordered 

Repair 


Count  Limit 


17297  Infinite 
2838  Infinite 
14428  Infinite 


Rim  Time:  1  hr(s)  53  min(s)  and  30  sec(s) 
Simulation  run  complete. 
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